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PREDiCfmG  THE  PERFORMANCE  OF  BAND  7 
COMMUNlCAf  ION  SYSTEMS 
USING  ELECTRONIC  COMPUTERS 


by 

Bonald  L^  Lucas  and  George  W,  Haydon 

ABSTRACT 

Radio  system  parameters  are  combined  with 
geophysical  and  ionospheric  cbaracteristics  to  pre¬ 
dict  the  performance  of  high  frequency  sky-wave 
Gommunication  circuits  through  the  use  of  electronic 
computers.  A  program  is  presented  to  compute 
Maximum  Usable  Frequencies »  Optimum  Traffic 
Frequencies,  Lowest  Useful  Frequencies,  probable 
mode  of  propagation,  angle  of  arrival,  circuit 
reliability,  system  loss,  available  signal-to-noise 
and  field  strength.  Numerical  representation  is 
used  for  all  parameters  not  expressed  in  closed 
mathematical  form. 

I,  INTRODUCTION 

Long  distance  high  frequency  radio  systems  are  subject  to 
marked  variations  in  performance,  most  of  which  are  directly  related 
to  changes  in  the  ionosphere.  Changes  in  the  ionosphere  result  in 
variations  in  the  maximiim  frequency  which  will  be  returned  back  to 
earth,  in  the  strength  of  the  radio  waves  due  to  increased  absorption, 
in  optimum  vertical  angles  of  wave  arrival  and  departure,  and  in  the 
background  noise  level  as  atmospheric  radio  noise  is  propagated  from 
great  distances.  Effective  operation  of  long  distance  high  frequency 
radio  systems  increases  in  proportion  to  the  ability  to  predict  variations 
in  circmt  performance  and  thereby  optimize  frequency  selection,  antenna 


ckoice  and  other  circuit  parameters  to  capitalize  on  anticipated  ionos¬ 
pheric  conditions.  It  is  the  purpose  of  this  report  to  present  a  current 
Computer  routine  which  has  been  developed  to  combine  the  more  pre¬ 
dictable  ionospheric  characteristics  with  circiut  parameters  to  calculate, 
the  expected  performance  of  high  frequency  sky-wave  radio  systems. 
Basic  ionospheric  characteristics  which  are  considered  predictable 
enough  to  be  useful  are: 

(1)  The  monthly  median  of  the  ordinary  ray  vertical  incident 
critical  frequencies  of  the  F2  layer  (foF2). 

(2)  The  monthly  median  of  the  factors  relating  the  vertical 
incident  critical  frequency  to  Maximum  Usable  Frequency  on  a  3000  km 
path  (M-3000  factor). 

(3)  The  monthly  median  minimum  virtual  height  of  the  F2  layer. 

The  following  ionospheric  characteristics  are  derived  in  the 

prediction  process: 

(1)  The  F -layer  critical  frequency  as  a  function  of  solar  activity 
and  angle  of  the  sun. 

(2)  Ionospheric  absorption  as  a  function  of  operating  frequency, 
gyro  frequency,  solar  activity,  angle  of  the  sun,  and  vertical  angle  of 
departure . 

(3)  Factors  relating  the  E-layer  critical  frequencies  at  different 
departure  angles. 

(4)  Factors  relating  F-layer  critical  frequencies  at  different 
departure  angles. 

The  following  geophysical  data  are  utilized: 

(1)  E-region  gyro  frequency. 

(2)  F-region  gyro  frequency. 

(3)  Geomagnetic  latitudes, 

(4)  World-wide  distribution  of  land  and  sea. 
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(5)  World-wide  distribution  of  atmospberic  noise  including 
seasonal,  diurnal  and  frequency  variations. 

(6)  Cosmic  noise  including  frequency  dependence. 

(7)  Man-made  noise  in  the  receiving  area  including  frequency 
dependence. 

The  following  circuit  data  are  needed: 

(1)  Transmitter  location. 

(2)  Receiver  location. 

(3)  Transmitter  power. 

(4)  Transmitting  antei^a  (physical  dimensions  and  orientation). 

(5)  Receiving  antetma  (physical  dimensions  and  orientation). 

(6)  Type  of  traffic  (modulation,  type  of  intelligence,  speed  of 
transmission,  required  quality). 

The  above  data  are  combined  to  determine  the  following  perfor¬ 
mance  characteristics  of  the  circuit: 

(1)  The  Maximum  Usable  Frequency  (MUF). 

(2)  The  Optimum  Traffic  Frequency  (FQT). 

(3)  System  loss  (ratio  of  power  received  to  power  transmitted). 

(4)  Field  strength  at  the  receiver. 

(5)  Monthly  median  of  hourly  median  available  signal -to -noise 

ratio. 

(6)  Circuit  reliability  (per  cent  of  days  the  hourly  median  signal 
to-noise  ratio  can  be  expected  to  equal  or  exceed  a  given  value). 

(7)  The  Lowest  Useful  Frequency  (LUF),  the  frequency  below 
which  the  circuit  reliability  is  expected  to  be  less  than  90%. 

In  addition  to  the  above  circuit  performance  characteristics,  the 
following  other  circuit  characteristics  are  determined: 

(1)  The  short  great  circle  distance. 

(2)  Bearing  of  the  transmitter  at  the  receiver. 


-  4  - 


(3)  Bearing  of  the  receiver  at  the  transmitter » 

(4)  Propagation  path  having  the  lowest  system  lossi  e.g.  *  2  hops 
via  F  layer  -  ZF* 

(I)  Vertical  reception  angle  associated  with  the  best  propagation 

path. 

The  Solution  of  the  problem  is  divided  into  two  parts:  (1)  an 
estimation  of  the  available  signal,  and  (2)  an  estimation  of  i^e  required 
signal.  The  available  signal  depends  upon: 

(1)  The  transmitter  power. 

(2)  The  transmitter  antetma  gain. 

(3)  Loss  due  to  spreading  of  the  radio  energy  as  it  propagates 
to  greater  distance. 

(4)  Losses  in  the  ionosphere  due  to  absorption. 

(5)  Losses  at  each  ground  reflection. 

(6)  The  gain  of  the  receiving  antenna. 

The  required  signal  depends  upon: 

(1)  The  atmospheric  noise  at  the  receiver  location^ 

(2)  The  man»made  noise  at  the  receiver  location. 

(3)  Cosmic  noise. 

(4)  The  required  signal  ■•to  moise  ratio  depending  Upon  type  and 
quality  of  service  desired. 

The  most  basic  calculation  in  the  estimation  of  available  signal 
power  is  the  system  loss  calculation.  The  ionospheric,  geophysical 
and  circuit  characteristics  are  combined  to  estimate  a  quasi  ••minimum 
system  loss  (the  minimum  loss  expected  on  any  day  of  the  month  at  the 
hour  for  which  the  calculations  are  made).  The  distribution  of  losses 
above  the  quasi  ••minimum  have  been  empirically  determined  and  are  a 
f\mction  of  geographic  location  and  time  of  day.  These  statistical 
distributions  are  used  to  account  for  the  numerous  other  factors  wMch 
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contribute  to  variations  between  predicted  Signal  levels  and  tbose  ob¬ 
served  on  a  given  day.  Tbese  factors  include  polarization  misniatck 
between  tbe  signal  and  tke  anteima,  focusing  by  tbe  ionosphere ,  defo- 
cusing  by  the  earthi  variations  between  theoretical  and  actual  antenna 
performance  and  day  to  day  variations  in  ionospheric  layer  heights^ 
ionospheric  absorption  and  critical  frequencies. 

Basic  F2-layer  data  for  long  term  predictiohs  are  numerical 
maps  of  the  F2  region  as  determined  from  world-wide  observations 
using  19S4  as  typical  of  a  low  solar  activity  period  and  1958  as  typical 
of  a  high  solar  activity  period.  In  addition  to  1954  and  1958,  F2  data 
from  other  periods  of  high  and  low  solar  activity  have  been  incorporated 
in  the  March,  June,  September  and  December  predictions,  since  these 
seasonal  extremes  are  often  used  in  long  term  planning. 

Numerical  maps  of  the  F2  layer  for  monthly  predictions  will  be 
a  part  of  the  "ionospheric  Predictions"  CRPL  Series  D  issued  by  the 
Central  Radio  Propagation  Laboratory  starting  with  the  January  1963 
issue. 

Sporadic  E  and  F  1-layer  predictions  are  considered  in  this  report, 
only  insofar  as  propagation  by  these  modes  enter  into  the  empirically 
determined  signal  distributions. 

II.  PROGRAM  DESCRIPTION 

The  computer  routine  is  based  on  established  manual  methods 
and  assiunes  a  working  knowledge  of  these  methods  [  Laitinen,  1949], 
[Haydon,  1962],  [NBS  Circular  No.  462],  along  with  a  familiarity  of 
Fortran  II  computer  language  [McCracken,  1962]. 

The  computer  method  described  is  designed  for  any  (32K  IBM- 
7090  class)  computer  with  the  necessary  compiler  and  tape  units.  The 
Fortran  II  computer  language  assures  flexibility  and  ease  of  modification 


as  better  knowledge  of  the  various  parameters  beGomes  available.  The 
program  provides  the  communicator  with  a  completely  general  method 
for  predicting  the  performance  of  any  sky-wave  high  frequency  radio 
system.  The  system  performance  calculations  may  be  performed  for 
any  month  of  the  year,  and  any  degree  of  solar  activity  for  any  hour  of 
the  day.  Numerical  mapping  is  included  for  those  parameters  not 
evaluated  directly  by  the  closed  mathematical  expressions  ( Section  IX). 

The  program  is  designed  to  be  used  with  either  a  "long-term" 
prediction  data  tape  which  utilizes  coefficients  foF2  and  M-BOOO  factors 
representing  periods  for  high  and  low  solar  extremes,  or  a  "short-term 
data  tape  which  utilizes  monthly  coefficients  issued  by  CRPL  approx¬ 
imately  three  months  in  advance. 

MUF-FOT 

The  Maximum  Usable  Frequency  and  Optimum  Traffic  Frequency 
are  based  on  F2-region  numerical  mapping  of  ionospheric  character¬ 
istics  recently  developed  by  CRPL  [Jones  and  Oallet,  1961],  and  a 
semi-empirical  relationship  between  the  sun' s  zenith  angle  and  solar 
activity  for  the  regular  E  region  [Haydon,  1962].  Conventional  pre¬ 
diction  methods  [NBS  Circular  No.  462]  using  control  points  along  the 
great  circle  2000  kilometers  from  each  terminal  for  the  F2  layer,  and 
control  points  1000  kilometers  from  each  terminal  for  the  regular  £ 
layer  are  used  for  paths  greater  than  4000  kilometers,  The  midpoint 
of  the  path  is  considered  as  the  control  point  for  both  £  and  F  layer  on 
paths  less  than  4000  km. 

Since  the  basic  data  for  the  predictions  are  monthly  median  values 
predictions  for  specific  days  within  the  month  are  not  available. 
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2i  Sy8terft_LQM 

1.  System  loss  computations  are  based  on  rayapatb  tbeory  witb 
patb  limitations  based  on  ray-patk  geometry  between  a  spberical  eartla 
and  a  Goncentric  ionospbere.  Limitations  fall  into  two  categories, 
arrival  angle  and  penetration  frequenGy,  and  tbe  probable  paths  for  a 
given  circuit  at  a  given  hour  and  frequenGy  are  chosen  on  the  basis  of 
these  limitations. 

F2 -region  propagation  is  considered  unlikely  in  these  computations 
if  the  angle  at  the  F2  region  is  too  sharp  for  F^layer  support  or  if  the 
penetration  frequency  of  the  E  region  rises  to  a  value  which  will  support 
the  operating  frequency,  in  this  last  case,  £  region  transmission  is 
geometrically  possible.  Multiple  hop  transmission  for  a  given  circuit 
is  possible  solely  by  the  E  region  or  the  F2  region,  but  the  computer 
routine  also  includes  the  case  where  some  hops  of  a  given  path  are 
supported  by  the  £  region  and  other  hops  supported  by  the  F2  region. 

Three  F2-region  paths,  two  £-region  paths,  and  one  £<i>F2-regiQn 
path  are  inspected  for  a  given  circuit  at  a  given  hour  for  a  given  Ire# 
quency.  The  path  that  is  geometrically  possible  with  the  least  theo¬ 
retical  loss  is  assumed  to  be  the  most  useful  mode  with  a  loss  typical 
for  the  circuit.  The  total  field  is  not  included  in  this  routine. 

2.  Layer  heights  used  for  the  F2  region  are  minimum  vertical 
heights  averaged  for  a  given  month.  £ -layer  height  is  assumed  to  be 
110  km.  The  average  height  of  the  absorbing  region  is  assumed  to  be 
100  km. 

3.  Ionospheric  layer  heights  at  the  reflection  points  are  averaged 
in  calculation  of  the  vertical  radiation  angle  (4)  and  the  angle  of  incidence 
at  the  ionosphere  (<!>')«  Unequal  take-off  angles  in  multi-bop  cases  are 
not  included. 
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4.  Absoi^ptidn  iii  the  lower  regions  is  represented  by  a  senii* 
enapirical  relationship  involving  the  angle  of  incidence  at  the  idnospherei 
the  gyro  frequency  in  tiie  absorbing  region,  the  zenith  angle  of  the  sun, 
Solar  activity,  and  the  operating  frequency. 

5.  Estimation  of  the  ground  reflection  losses  is  evaluated  for 
either  land  or  sea  reflections.  Shore  line  and  mixed  paths  are  not 
included  and  random  polarization  of  the  downcoming  waves  is  assumed. 

6.  Convergence  and  divergence  of  the  ionospheric  waves  are 
inGluded  in  the  routine  only  in  the  empirical  determined  system  loss 
distributions . 

7.  Theoretical  vertical  and  horizontal  reflection  coefficients 
are  evaluated  for  the  ground  losses  and  the  antenna  patterns. 

8.  Receiving  antenna  response  (gain)  is  assumed  to  be  equal  to 
the  antenna  gain  of  an  identical  antenna  used  as  a  transmitting  antenna. 
All  gains  are  relative  to  an  isotropic  in  free  space. 

The  complete  computer  routine  contains  a  subroutine  for  three 
general  types  of  antennas. 

a.  Terminated  rhombics  at  any  height  above  the  terrain,  any 
leg  length  and  tilt  angle.  Poor  earth  (O'  =  .  001  ndios/  meter;  c  =  4)  is 
assumed  as  the  terrain,  but  any  ground  conductivity  and  dielectric 
constant  could  be  used.  Off -azimuth  radiation  may  be  computed  for 
the  rhombics. 

b.  Horizontal  half-^wave  dipoles  for  any  height  above  finite  earth, 

c.  Vertical  ship  antennas  for  constant  length  elements  or  for 
multiples  of  any  wave  lengths  not  less  than  ,  OZ  Practical  efficiencies 
of  grounded  verticals  less  than  one'^’quarteF  wave  length  are  computed 
and  included  in  the  antenna  gain. 


! 

i 
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9^  Signal  diatribution  is  assumed  to  vary  with  geomagnetiG 
latitude  and  length  of  the  eirGuit  and  quasi -minimum  losses  are 
adjusted  to  median  losses  in  terms  of  these  empiriGal  distributions 
[Haydont  1962]. 

10.  Off-great-GirGle  transmissions  are  not  inGluded  in  these 
GalGulations. 


3.  Field  Strength  Qgtion 

The  field  strenglh.  option  is  Gomputed  similar  to  the  system  losSi 
then  mathematiGally  related  to  the  system  loss  (Section  IX). 

4.  Radio  Noise 

The  CCiR  World-Wide  maps  of  atmospheric  noise  are  used  as 
representative  of  the  distribution  of  noise  [CCIR  Report  No.  65]. 

Man-made  noise  estimates  are  based  upon  the  type  of  receiving 
arSa  iHaydon,  1962].  Measured  values  of  man-made  noise  are  also 
acceptable. 

Galactic  noise  at  the  receiver  site  [CCIR  Report  No.  65]  is 
estimated  by  .a  least  squares  polynomial  in  frequency.  The  highest  of 
the  three  noise  types  is  assumed  to  predominate  and  the  total  noise 
field  is  not  estimated.  All  noise  values  are  expressed  in  db  relative 
to  1  watt  for  a  one -cycle  bandwidth.  The  1  Me/  s  atmospheric  noise 
above  KTB  is  adjusted  to  a  1  c  /  s  band  relative  to  one  watt  by  least 
squares  polynomial  representation  of  the  seasonal  maps  [Lucas  and 
Harper,  1962]. 

5.  Available  Signal « T  o  -No is e 

The  monthly  median  of  the  hourly  median  available  signal  is 
compared  with  the  average  noise  in  a  four<^hour  time  block  at  the 
receiver  site.  This  gives  an  estimate  of  the  monthly  median  of  hourly 
median  total  signal  po^er  relative  to  noise  in  a  one -cycle  band. 
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6.  Required  Signa.l-To*Ndi#e 

The  required  signal-to-noise  depende  upon  the  service  required^ 
e.  gi ,  error  rate  in  teletype^  intelligibility  of  voice  *  etc^  A  basic 
signal vto “noise  requireinent  based  upon  a  steady  static  signal  and 
random  noise  adjusted  for  signal  fading  degradation,  diversity  improve <• 
ment,  etCi  ^  determines  the  required  signal  “to  *noise. 

7i  Reliability 

Reliability  is  estimated  at  discrete  frequency  bands  by  calculating 
the  ratio  of  the  available  monthly  median  signal  “to  “noise  of  Paragraph 
S  to  the  required  hourly  median  signal“to“noise  of  Paragraph  6.  Using 
the  above  parameter  and  the  operating  frequency^  theoretical  reliability 
is  obtained  for  a  given  circuit  of  a  given  length  and  geomagnetic  latitude 
by  use  of  statistical  two  “dimensional  maps. 

8*  Lowe  St  U  seful  Frequency 

The  Lowest  Useful  Frequency  is  calculated  by  stepping  through 
the  HF  band  until  a  reliability  of  90%  or  more  is  obtained,  then  stepping 
down  in  small  increments  of  frequency  until  a  frequency  is  found  that 
has  (within  the  tolerance  of  the  frequency  increment)  90%  reliability. 
This  frequency  is  the  classical  90%  LUF. 

in.  PROGRAM  USE 

The  IBM“7090  Fortran  II  listing  is  shown  in  Section  VIU.  The 
cards  corresponding  to  these  listings  with  the  associated  control  cards 
necessary  for  the  specific  monitor  being  used  are  all  that  is  necessary 
to  compile  the  program,  A  pictorial  view  of  a  sunple  deck  ready  to 
run  under  the  1BM“7090  monitor  is  shown  in  Figure  17.  System  sub“ 
routmes  used  in  the  program  are  as  follows: 
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Subroutine 

Call 

Argument 

1 .  Sine 

SINF  (x) 

(radians) 

2.  Arc  sine 

ASINF  (x) 

(-  1  S  X  «  1) 

3.  Cosine 

COSF  (X) 

(radians) 

4.  Arc  cosine 

ACOSF  (x) 

(*  1  ^  X  «  1 ) 

5.  Absolute  value 

ABSF  (x) 

(floating ) 

6.  Absolute  value 

XABSF  (x) 

(fixed) 

7.  Arc  tangent 

ATANF  (x) 

(i.  00  ^  X  ^00 

8.  Square  root 

SQRTF  (x) 

(0  ^  X  ^00) 

9.  Transfer  of  sign 

SIGNF  (x,y) 

(arg  L  arg  2) 

10.  Logarithm  (Base  10) 

LOCI  OF  (x) 

(0  <  X  ^  ec  ) 

11.  Natural  logaritlun 

LOGF  (X) 

(  0  <  X  ^  ee ) 

12.  hiteger  function 

INTF  (x) 

(.  00  ^  X  ^  oe 

The  entire  prograni  including  common  variables  and  system 
Subroutines  occupies  approximately  24000  decimal  locatipnSi 

Boutines  included  in  the  program  with  a  brief  description  of  their 
function  are: 

MUFLUF  --  (2550  decimal  locations) 

This  is  the  main  program.  It  calculates  the  Maximum  Usable 
Frequencies,  Optimum  Traffic  Frequencies,  and  other  parameters 
needed  not  only  in  the  calculation  of  these  critical  frequencies  but  also 
in  one  or  more  subroutines. 

LUFFY  (3857  decimal  locations) 

This  subroutine  is  larger  and  more  complex  than  MUFLUF,  but 
is  a  subprogram  because  the  option  is  asked  for  in  MUFLUF.  This 
routine  calculates  the  systems  performance  of  a  given  circuit,  i,  e, , 
system  loss,  reliability*  signal«to^noise  and  field  strength.  The  subn 
programs  that  follow  are  used  by  LUFFY  to  perform  a  specific  task 
needed  in  the  system  performance  calculations. 
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POLY  --  (107  decimal  locations) 

Tkis  snbf outine  is  called  upon  to  evaluate  all  parameters  that 
are  represented  by  Nth  degree  polynomials  of  the  power  series  in  X. 


Example: 

Y  =  (As  s  X  ^  +  As  ,  X 

V  0,0  0,1 

1  .  .A  ft  ^  ^  0 

4^  ♦  •  •  +  A.  3C  j  Z 

O^n  y 

• 

...+ As  x  "^  Z  ^ 
l,n  y 

• 

+  i  s  • 

^(a  x°  +  A 

\  m,0  m,l  m,n  y 

LOSS  -«  (214  decimal  locations) 

This  subroutine  evaluates  the  ground  reflection  losses  as  a 
fmiction  of  the  reflection  angle  and  operating  frequency  over  finite 
earth.  Random  polarization  is  assumed  in  the  calculations  as  the 
vertical  and  horizontal  reflection  Coefficients  are  assumed  to  contri¬ 
bute  equally  in  the  polarization  effects  of  the  loss. 

This  routine  is  unique  in  that  complex  arithmetic  is  absent  in 
the  calculations  of  these  reflection  losses  [  Phillips.  19^1]. 

GAIN  --  (684  decimal  locations) 

This  routine  calculates  theoretical  power  gains  of  horizontal 
rhombics  in  three  dimensions .  horizontal  dipole s.  and  practical 
grounded  vertical  radiators,  all  over  finite  earth.  Practical  lower 
limits  of  gain  have  been  set  as  -10  db  relative  to  an  isotropic  in  free 
space. 

Cl  9-  (411  decimal  locations) 

This  fvmction  subroutine  evaluates  the  sine  and  cosine  integral 
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X 

Si  (x)  =  ^  \  sin  V  dv 
0  V 

X 

Gi  (x)  =  Q  C  dv 

*  V 

wken  called  for  by  the  loss  and  gain  routines. 

NOISY  --  (215  decitnal  locations) 

This  routine  calculates  the  world-wide  distribution  of  atmospheric 
noise  and  major  land  bodies  which  are  represented  by  Fourier 
coefficients. 

The  evaluation  is  made  at  a  given  longitude  and  latitude  of  the 
receiving  location  or  ground  reflection  points 
CURVY  --  (1462  decimal  locations) 

This  routine,  although  not  pertinent  to  the  solution  of  any  problem 
encountered  in  system  performances,  produces  a  semi-line  graph 
representation  of  the  diurnal  variation  of  either  or  both  the  Optimum 
Traffic  Frequency  and  the  Lowest  Useful  Frequency. 

VREFCO  --  (230  decimal  locations) 

This  routine  calculates  the  ground  reflection  coefficient  (K^  )  for 
vertically  polarized  waves. 

NOTE; 

The  ground  reflection  coefficient  (K  )  for  horizontally  polarized 

w 

waves  is  calculated  within  subroutine  GAIN. 

IV.  DESCRIPTION  OF  DATA  TAPES 
1.  "Long-term"  Data  Tape. 

The  binary  data  tape  for  the  storage  of  coefficients  representing 
those  variables  not  expressed  in  closed  form  has  a  packing  density  of 
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556  frames  per  inch.  Twelve  logiGal  records  -  one  for  each  month  « 
are  used.  Each  logical  record  is  composed  of  physical  redords  256 
Words  in  length.  Approximately  800  feet  of  tape  are  used. 

Each  logical  record  contains  six  one  >dimensionar  arrays «  three 
three-dimensional  arrays^  one  two-dimensional  arrya  in  the  following 
order: 

IL  (4),  JL  (4).  KL  (4),  LK  (4),  JAL  (4),  Q  (20,60,4)* 

A  (10,  7,  14),  P  (29.16,6),  ABP  (2,  6) 

The  variables  IL,  JL,  KL,  LK,  and  JAL  are  limits  to  which  the 
Fourier  generation  of  the  foF2  and  M-3000  factors  are  parried,  e^  gi , 
|L  (1)  is  associated  with  the  two  dimensional  set  of  coefficients 
Q  (M.N.l). 

The  "Q"  arrays  [Lucas,  1961]  contains  the  coefficients  used  in 
the  generation  of  the  foF2  and  M-3000  factors  in  the  following  order 
for  each  month: 

1.  Q  (M,N,  1)  foF2  -  low  Solar  activity 

2.  0  (M,  N,  2)  foF2  -  high  solar  activity 

3.  Q  (M,  N,  3)  M-3000  -  low  solar  activity 

4.  Q  (M,  N,  4)  M-3000  -  high  solar  activity 

The  "A"  array  contains  the  coefficients  needed  for  other  para¬ 
meters  not  represented  in  closed  form* 

A  (M,N,  1)  gyro -frequency  distribution  [  Lucas,  1961} 

A(M,  N,2)  H'f2  layer  height  charts  [Haydpn,  1962] 

A  (M, N,3)  blank  to  leave  space  for  numerical  representation 
of  measured  patterns  of  antenna.s 
A  (M,  N,4)  nighttime  frequency  dependence  of  atmospheric 
noise  [Lucas  and  Harper,  1962], 

A  (M,  N,  5)  daytime  frequency  dependence  of  atmospheric  noise 
A  (M,  N,  6)  man-made  noise  level  in  industrial  area  [Haydon, 

I  962] 
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A  (M^N,  7)  man-made  noise  level  in  residential  area 
A  (M^N,  8)  man-made  noise  level  in  rural  area 
A  (M,N,9)  man-made  noise  level  in  remote  luipopulous  area 
A  (M^N,  10)  daytime  distribution  of  circuit  reliability  in 
temperate  regions  [Haydoni  1962] 

A  (M,Ni  11)  nigbttime  distribution  of  Circuit  reliability  in 
temperate  regions 

A  (MiNi  12)  distribution  of  circuit  reliability  in  polar  region 
A  (KiiNi  13)  distribution  of  short  circuit  reliability  in 
auroral  region 

A  (M,N,  14)  distribution  of  long  Circuit  reliability  in 
auroral  region 

The  "P"  array  of  each  logical  record  contains  five  maps  of  the 
world«»>vide  distribution  of  atmospheric  noise  [Lucas  and  Harper,  1962] 
for  the  month  being  Gonsidered.  The  major  land  bodies  of  the  world 
are  also  represented  in  this  array. 

P  (M,N,  1)  2000-0400  LMT 

P  (M,N,2)  1600-2000  LMT 

P(M,N,3)  1200-1600  LMT 

P(M,N,4)  6800-1200  LMT 

P(M.N,  5)  0400-0800  LMT 

P  (M»N,6)  major  land  bodies  of  the  world 
The  array  "ABP"  contains  the  convergence  factors  for  the  Fourier 
generation  of  the  "P"  matrix  [Lucas  and  Harper,  1962]. 

The  elements  of  the  matrices  that  are  not  sensitive  to  month  are 
identical  for  all  logical  records. 

2,  "Short-term"  Data  Tape. 

(a)  Description  of  "Short-term"  Data  Tape 
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The  "short-term"  data  tape  is  identical  to  the  "long-term"  data 
tane  v.ith  the  exception  of  the  £oF2  and  the  M-3000  coefficients.  On 
the  "short-term"  tape  each  particular  month's  coefficients  are  stored 
in  positions  used  by  the  coefficients  representing  high  and  low  solar 
activity  extremes  on  the  "long  term"  tape.  The  solar  activity  index 
associated  with  the  "Short-term"  Coefficients  must  be  used  when  pre¬ 
dictions  are  made  using  these  Coefficients.  The  routine  v.hich  follows 
is  designed  to  generate  and  update  a  "short-term"  data  tape  as  the 
coefficients  become  available.  It  is  designed  to  use  the  "long-term" 
or  "short-term"  tape  (Logical  »2  and  Logical  tape  *3)which  will  become 
the  new  "short-term"  data  tape.  Precede  the  data  deck  w-ith  one  card 
which  contains  the  first  and  last  months  for  which  the  "short-term" 
tane  is  to  be  uDdated,  A.  minus  1  Card  punched  in  Columns  55-36  must 
follow'  each  individual  £oF2  or  M-3000  coefficient  deck.  A  sample  deck 
for  one  month' s  data  is: 


• ::ra  s  X  L  r  ^  :s r  »  r  ansi ‘t  »:cM6  *  r  i: 5.U  !» F  u Si- 1:  c -i  !*  vi « *:*  *:;“■»  >  ^ 

1 1 1  n  1 1 11 1 1 1 1  H  n  1  n  !  M  n  M  i  n  nj 


1 1 11 !  1 1  n  n  i  1 1  m  M  n  n  n  n  11 1 1  n  n  n  1  n  I  n  1 
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(b)  Program  to  Use  GRPL  Series  D  Predicted  Goefficients 
to  Generate  and  Update  "Sbort-term''  Tape. 

c  always  use  most  recent  tape  for  updating  SM0RT*TERM  tape 

OlMENSrON  Q(2d*60«4)»!L(4)«JL(4)«iCL(4)«LK(4)»JAL(4)*Aa6*7fU)fPl( 

i2^»l6t6)*AiP(2i6) 

1  FORMAT  (l9X*2l4*4X»2t4tiaXt2lStEl7*7) 

2  FORMAT  (5iX*2I5«E17«7) 

101  F0RMAT(2I2) 

RIAO  INPUT  TAPE  5tl01tM0NTMl»M0NTH2 

REWIND  2 

rewind  3 

MONTH3*MONTH2-»-i 

LOCK  ■  M0NTM1«1 

IF (LOCK)  4t55*4 

4  DO  3  II»l*LOCK 

READ  TAPE  2* tLtJL*KL*LK*JALfO»A»Pl»ABP 

3  WRITETAPE  3* lLtJL*KL«LK»JAL»0»A»Pl*ABP 
B9  DO  lOQ  II«MONTH1»MONTM2 

READ  TAPE  2*  !LfULtKL»LK»JAL>Q»A»Pl*ABP 
10«l 

5  READ  INPUT  TAPE  5»l»I»J»KfL»M>N»P 
rL(IO)«I'»'l 

JL(|0)>J<fl 
KL(lO)«K4l 
LK(I0)«L<fl 
JAL(  |0)92«J<i>1 
20  M»M4.1 
M»N+1 

Q(MfNtIO)«P 

READ  INPUT  TAPE  9t2»MrM»P 
|F(M)  40»20»20 
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40  1 0»  10^*2 

If (4-10)  6*5*5 

6  DO  7  I0>a*4*2 
tL(lO)slL( 10-1) 

JL(iO)»Jt.(  10-1) 

KU(I0)*XL« I0“1) 

LK(  iO)»LIC(  lO-l) 

JAL(I0)*JAL( 10-1) 

DO  7  M«1*2Q 

DO  7  N«l*60 

7  Q(M*N*10)>0(M*N*!0-1) 

100  WRlTEtAPE  3*  IL* JL*KL.LK* JAL*0* A»P1 *ABP 

1F(M0NTH3-12)  301*30lf40l 
301  DO  201  K«MONtH3«12 

READ  TAPE  2*IL*JL*KU*L<*JAU*Q*A*PI*ABP 
201  WRITETAPE  3*iL*JL*Ki*LK*JAL«0*A*Pl*ABP 
401  REMIND  2 
REWIND  3 
CALL  EXIT 
END 


Minor  changes  in  the  Fortran  program  found  above  wovild  allow 
the  updating  of  the  "long-term"  tape  as  better  coefficients  become 
available. 

E^reme  caution  should  be  exercised  when  changes  are  being 
made  to  the  data  tape.  Fortran  statements  identical  to  the  read  state 
ments  in  the  main  program  "MUFLUF"  should  be  used  to  re  ••write  the 
tape. 

Inspection  and  checking  of  coefficients  contained  on  the  tape  may 
be  accomplished  by  merely  dumping  the  tape  off-line.  It  is  suggested, 
however,  that  for  easier  inspection  it  be  read  by  the  "read  tape"  state¬ 
ments  contained  in  "MUFLUF", 


Data  tapes  for  the  IBM^7O90  containing  all  the  logical  records 
may  be  obtained  from  the  Radio  Systems  Division,  National  Bureau  of 
Standards,  Boulder,  Colorado. 

V.  PROGRAM  OPTIONS 

Options  inGluded  in  the  program  aid  in  the  solution  of  a  wide  range 
of  systems  problems.  Options  designed  to  be  printed  on  standard 
9  by  11  paper  are: 

1.  MUF  and  FOT  (Figure  1) 

The  Maximum  Usable  Frequency  (MUF)  and  Optimum  Traffic 
Frequency  (FOT)  are  printed  three  circuit  months  per  page. 

2.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE,  AND  RELIABILITY 
(Figure  2) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  circuit 
reliability  are  printed  one  circuit  month  per  page.  MUF*FOT'  s  are 
calculated  for  each  hour  of  GMT.  Mode,  angle  and  reliability  are 
calculated  only  for  even  hours  of  GMT. 

3.  FOT  and  LUF  (Figures  3  and  4) 

The  Optimum  Traffic  Frequency  and  the  Lowest  Useful  High 
Frequency  are  printed  three  circuit  months  per  page.  The  Optunum 
Traffic  Frequency  is  always  given  each  hour  of  GMT  while  the  cORmnun. 
icator  has  a  choice  of  each  hour  or  every  even  hour  of  GMT  for  the 
Lowest  Useful  High  Frequency. 

4.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE  AND  SYSTEM  LOSS 
(Figure  5) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  monthly 
median  of  hourly  median  system  loss  are  printed  one  circuit  month  per 
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page.  MUF  and  FOT'  s  are  caiLculated  for  all  hoinrs  of  GMT  While  modei 
angles  and  system  loss  are  calculated  for  even  hours  of  GMT. 

5.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE  AND  FIELD  STRENGTH 
(Figure  6) 

The  Maximum  Usable  Frequency  and  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  monthly 
median  of  the  hourly  median  field  strength  are  printed  one  circuit  month 
per  page.  MUF  and  FOT*  s  are  calculated  for  all  hours  of  GMT  while 
mode,  angle,  and  field  strength  are  calculated  for  even  hours  of  GMT. 

6.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE  AND  SlGNAL-TO* 

NOISE  (Figure  7) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
probable  mode  of  propagation,  theoretical  angle  of  arrival  and  monthly 
median  of  the  hourly  signal i-to -noise  at  the  receiving  antenna  terminals 
are  printed  one  circuit  month  per  page.  MUF  and  FOT*  s  are  calculated 
for  all  hours  of  GMT  while  mode,  angle,  and  signal-to-noise  are  cal¬ 
culated  for  even  hours  of  GMT. 

7.  MUF,  FOT  AND  reliability  (Figure  8) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency 
and  circuit  reliability  are  calculated  every  even  hour  of' GMT  and 
printed  two  circuit  months  per  page. 

8.  FOT  GRAPHS  (Figure  9) 

The  Optimum  Traffic  Frequencies  are  plotted  frequencies  vs. 

GMT  time  on  a  nonlinear  frequency  scale  by  the  on-line  printer. 

9.  FOT  AND  LUF  GRAPHS  (Figure  10) 

The  Optimum  Traffic  Frequencies  and  Lowest  Useful  High 
Frequencies  are  plotted  frequencies  vs.  GMT  time  on  a  nonlinear 
frequency  scale  by  the  on-line  printer. 
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The  above  options  may  utilize  one  of  two  input  data  tapes.  Long« 
term  predictions  are  available  from  a  data  tape  with  numericai  coeffici* 
ients  for  the  months  representing  both  high  and  low  solar  activity. 

These  extremes  are  used  to  predict  MUF  associated  with  the 
intermediate  sunspot  number. 

Short-term  predictidns  are  available  from  a  data  tape  with 
numerical  coefficients  for  approximately  three  months  in  advance.  No 
interpolation  is  made  when  using  this  data  tape. 

VI.  COMPUTER  INPUTS 


1  i  Circuit  Card  Format 


A  data  sheet  should  be  prepared  for  each  circuit  before  punching 
the  card.  The  number  in  the  left  hand  column  below  refers  to  infor¬ 


mation  on  the  data  sheet  (Figure  11).  The  fields  should  be  punched  as 
follows  (Figure  13)  with  numbers  right  justified. 

j^ta.Sheet  Card  Columns 


1. 

2. 

3, 

4, 


1-5 

6-11 

12-16 

17-22 


5. 

6. 

7. 

8. 

9.  (See  Note  #1) 
10.  (See  Note  #2) 
ll» 

12. 


28-30 

31-33 

34-36 

37-39 

40-41 

40-41 

42-44 

45-47 

48-50 


13. 
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14.  (See  Note  #1)  §U52 

15.  (See  Note  #2)  5U52 

,16.  (See  Note  #4)  53.55 

17.  (See  Note  #3)  53*55 

18.  56*58 

19.  Does  not  apply 

20.  (See  Note  #5)  59-62 

21.  (See  Note  #5)  63*66 

22.  71*72 
NOTES; 


1.  (a)  if  rhombic  is  circled,  enter  *1 

(b)  if  X/ 2  dipole  is  circled,  enter*3 

(c)  if  vertical  is  circled,  enter  *2 

(d)  if  vertical  in  multiples  of  wave  length 
is  given,  enter  in  height  cols.  *14  for 
X/ 4;  *12  for  X/2;  -1  for  X,  etc. 

2.  if  this  nunaber  is  given  instead  of  the 
information  which  appears  in  Note  1  (d) 
above,  enter  this  number. 

3.  (a)  If  check  in  Industrial,  enter  *1 

(b)  If  check  in  Residential,  enter  *2 

(c)  If  check  in  Rural,  enter  ^3 

(d)  If  check  in  Remote  Unpopulous,  enter  *4 

4.  If  number  is  given,  enter  this  number  ini 
lieu  of  #17  from  data  sheet. 

5.  Bearing  in  degrees  east  of  north. 
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2.  Met  jro  d  C  ar  d  F  or  m  at 

Tke  followiiig  is  a  Ust  of  methods  available <  The  numbers  on 
the  right  call  the  various  methods  when  placed  in  the  columns  shown. 


CardGol^^ns  Print* 

Galculatioas  to  be  Performed  2  4  6  8  Out 

MUF*FOT  listing  from  Gircuits  on  cards  110  0  Figure  1 

MljF*FOT  listing  from  Circuits  on  tape  10  0  0 

FOT  Curve  from  circuits  on  cards  110  1  Figure  9 

FOT  curve  from  circuits  on  tape  10  0  1 

Reliability *mode*angle  from  circuits  2  10  0  Figure  2 

on  cards 

Reliability-mode*angle  from  circuits  2  0  0  0 

on  tape 

FOT*LUF  listing  every  hour  from  3  10  0  Figure  3 

circuits  on  cards 

F0T*LUF  listing  every  hour  from  3  0  0  0 

circuits  on  tape 

FOT  listing  every  hour»  LUF  every  2  3  110  Figure  4 

hours  from  circuits  on  cards 

FOT  listing  every  hour,  LUF  every  2  3  0  1  0 

hours  from  circuits  on  tape 

« 

FOT*LUF  Curves  from  circuits  on  cards  3  10  1  Figure  10 

F0T*LUF  curves  from  circuits  on  tape  3  0  0  1 

System  loss  from  circuits  on  cards  4  10  0  Figure  5 

System  loss  from  circuits  on  tape  4  0  0  0 

Field  strength  (dbu)  irom  circuits  on  4  110  Figure  6 

cards 

Field  strength  (dbu)  from  circuits  on  tape  4  0  10 

Reliability*FOT  every  2  hours  from  5  10  0  Figure  8 

circuits  on  cards 

Re  liability  «FOT  every  2  hours  from  5  0  0  0 

circuits  on  tape 
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Available  S/  M  from  eireuits  on  ca^ds  6  10  0  Figure  7 

Available  S/  N  from  circuite  oa  tape 

3*  Assembling  the  Data  Decks 

^ly  one  data  tape  is  needed  with  these  routines  (Tape  #2).  A 
utility  tape  #3  is  alto  needed  wbicb  will  contain  the  circuits  from  a 
given  ruOi  If  a  large  number  of  circuits  are  to  be  run  repeatably,  it 
is  wise  to  save  the  tape  and  run  from  it  in  lieu  of  circuit  cards ^ 

To  assemble  a  data  decki  first  make  the  method  and  frequency 
complement  cards  for  l^e  desired  HF  predictions^  If  the  circuits  are 
to  be  read  from  cards,  the  circuit  cards  must  follow  the  frequency  card. 
Following  the  last  circuit  card.,  insert  a  "nines  card"  (the  number  9 
puncbed  in  columns  1  through  19),  to  indicate  the  end  of  the  circuit 
cards.  The  next  cards  will  indicate  the  months  and  solar  activity 
levels  for  which  the  calculations  are  to  be  made. 

The  months  January  through  December  are  indicated  by  niunbers 
1  through  12  placed  on  the  card  in  columns  23  and  24»  On  the  same 
card,  in  columns  25,  26  and  27,  punch  the  appropriate  sunspot  number. 
If  it  is  desired  to  terminate  the  calculations  of  these  circuits  under  the 
called  method,  a  card  with  ^1  punched  in  columns  23  and  24  should 
follow  the  "month  and  sunspot  number"  card.  However,  if  another 
method  is  to  be  run  for  the  same  circuits,  follow  the  -^1  card  by  method 
and  frequency  cards  that  allow  the  circuits  to  be  read  from  tape  (see 
method  card  format).  The  circuits  will  be  read  from  tape  #3  just 
prepared.  Following  the  new  frequency  and  method  cards,  place  the 
desired  month  and  sunspot  number  cards.  Notice  the  nine' s  card  has 
been  oimtted,  A  -I  card  mserted  after  the  month  and  sunspot  number 
card  will  terminate  the  method  and  a  blank  card  following  the  »l  card 
will  terminate  the  job. 
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The  data  deck  (to  he  rtm  from  cards)  is  composed  of  the  followmg 
cards  in  the  order  they  appear. 

1 .  Method  card 

2.  Frequency  complement  card 

3.  Group  of  circuit  cards 

4.  Nines  card 

5.  Group  of  month  and  svinspot  cards- 

6.  card 

(See  Figures  14  and  15) 

Any  number  of  different  methods  with  different  circuits  may  be 
set  up  sequentially  in  the  above  order.  When  it  is  desired  to  terminate 
the  job  I  place  a  blaidc  card  after  the  -1  card. 

A  data  deck  to  be  run  from  tape  is  composed  of  only: 

1.  Method  card 

2.  Frequency  complement  card 

3.  Group  of  month  and  sunspot  cards 

4.  *1  card 

(See  Figure  16) 

4.  Furpose  of  Cards  in  Data  Deck 

I 

(a)  Method  Card: 

The  method  card  shown  in  Figure  13  allows  the  communicator  a 
choice  of  computations  shown  in  the  example  print-outs  (Figures  1 
through  10)  with  the  option  of  having  data  for  the  circuits  on  either 
cards  or  magnetic  tape. 

(b)  Frequency  Complement  Card: 

The  frequency  card  allows  the  communicator  a  choice  of  10  fre¬ 
quencies  for  which  calculations  will  be  made,  e.g. ,  3.  1240,  4,  7856 
. .  30.0000.  tf  fewer  than  10  frequencies  are  wanted  merely  put 
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the  desired  ones  sequentially  on  card  from  the  begi^ing,  then  punch 
unused  fields  with  996000.  If  method  #3  is  to  be  run,  it  is  suggested 
that  a  frequency  card  identical  to  Figure  12  be  used. 

(c)  Circuit  Cards  (Figure  13): 

The  circuit  cards  are  designed  to  contain  all  necessary  circuit 
parameters  to  be  run  under  any  choice  of  the  computations  shown  in 
Figures  1  through  10. 

(d)  Nines  Card  (Figure  12): 

The  nines  card  is  a  circuit  card  which  is  not  computed.  The 
machine  senses  this  card  and  is  told  the  entire  list  of  circuits  has  been 
calculated  for  the  given  month. 

(e)  Month  and  Sunspot  Number  Card  (Figure  l2): 

The  month  card  contains  the  month  of  the  year  and  its  associated 
Sunspot  number  for  the  desired  circuit  calculations i 

(f)  Minus  One  Card  (Figure  12): 

The  minus  one  card  instructs  the  machine  that  all  computations 
are  complete  for  a  given  method. 

(g)  Blank  Card  (Figure  12): 

A  blank  Card  following  a  "minus  one  card"  instructs  the  machine 
that  all  computations  for  a  given  run  are  complete. 
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Figure  1.  Computer  Print  «Out  of  Circuit  MUF  and  FOT 
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Figure  8.  Alternative  Computer  Print^Qut  of  Circuit  MUF,  FOX  and 
Circuit  Reliability 
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TransmittinE  Antenna  Gain  if  Tyne  Other  Than  item  #9 


Receiving  Antenna  -  Height 

-  Leg  Length 

-  Tilt  Angle 

B  Type;  Rhombic,  \/  2  Dipole,  or  Vertical 
Receiving  Ahtetma  Gain  if  Type  other  Than_item.#  14 


Man-Made  Noise  at  Receiver  input 


Meters 

Meters 


Frequency 


Sandwidth 


DB  <  1  Watt 


17.  Type  of  Receiving  Area  if  #16  is  Unknown 


Industrial 


Residential 


Rural 


Remote  Unpopulous 


18,  Required  Hourly  Median  Signal -to -Noise  in  1  c/  s  Bandwidth 


19.  Detailed  Description  of  Serviceli  #18  is  Unknowir 

Type  of  Service f  S.g.  .  Radiotelephone  . . 

Modulation,  e.g.,  SSB  ^ 

Bandwidth  7  T  Tf' 

Nvunber  of  Channels  ~ 

Words  Per  Minute  Per  Channel  y  ~~ 

Maximum  Teletype  Error  Rate  '  ~  7, 

Minimum  Radiotelephone  intelligibility 
Diversity  Employed 

Any  Other  Description  of  Service  ^  ] 


20,  Enter  Bearing  of  Transmitting  Rhombic  Anterha  if  Off -Great- 

Circle  Pa^ 

21,  Enter  Bearing  of  Receiving  Rhombic  Antenna  if  Off-Great- 

Circle  Path 

22,  Two  Letter  Code  for  Receiver  Location 


Figure  1 1 
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Nines  Card 


w5i!!ttSir)i 

KUTtUtHt  -I 


NUMCH 

OiO'OODl 

4i|)  >  t  I 


f  .■‘ay^'.’^s .',  -■.<(’?- 
SiSf  ':Sii-  i'  j .  *  .!»-■  -•• 

FORTRAN  STATEMENT 
'too  «'d  d  d  a  dd  oTo  dB'tfTirnrrdiTiTnTrdiTiTdf  ddm"frd"f  dTd  cTTfoiTrc  odd  d 

•  •  U  II  nil  HI!  II II  tt  null  aill4B)illliS«lIBUMI>II>lll»»«  tIUWM  ■JSIUHKU  (leCcK  r.dill!  X 

.  . .  •  .  ■  . 


N  I  It  .  .1.1 1  ih  1 

5  0  0  0  c  Oddi 

h  K  c  ii ;/  7£  re 


Mdnth  and  Sunspot  Card 


lirAli'MENI 


NuMOkH 

\mn  bl 

ilj  1  4  i 


FORTRAN  STATEMENT 

c  0  0  0  0  u  Q'Q  0  (1 0  o"c  8  0  0  fi  0  fl  0  8  ?  0  8  0  8'fi  FFOTO  mTTtfym  FiTmTTim 


IwwvwuuuwuMWWVUUUVWWVUUUWWWWVVVVVWWVHVVWWVWUVWVVVVVVIdWUVHtfUVfVWWVVV 
jl  I  •  le  II  ’1  n  It  ti  i(  II  II  te:o  N  n  21N  )■  »N]|  n  »1I  tf  ]}  M  25  »  Jr3»]l4l4l  0«1  4t«6  4t.4;  #<;  !i  SI  -J  U  M  S7  SI  M  UU  W  KillC  CttttTC  )'  77 


ItyENTjf  r,d."10M 

figtrreF 

Its  M  n  R  V  n  n  loi 


dlo 


Month  =  3;  Sunspot  Number  14 


Minus  One  Card 


C^-foa 

5r&-TiLll£HT 

NUlMCA 

diinraj 

d 

111  9  4  9 

< 

«1 


■  > 


FORTRAN  STATEMENT 

fod  d'd  ff  d  0  odd  ii  i!  3  0  d  (!dTTr'ddTifinr'(rrdinnnriiTdd()  oTontTo  oii'3"iiii  ootb  cin'o  c  oTo  odoiiDTf  s  {ToI 

I  I  I  nil  tT  11 1(  li  It  IT  li  '.in vi;;  Tijtnntiaflaiis  uMKitiiaaMSMiMMad  tr-tr  st •:  y.-.t  xn  sec  m  ;?  u ti  c  ubt  nstt.  wn  /<  t  »  tstc:*  nnw 


FtequenGy  Complement  Card 

Moqs  ysoooo '  7(^  ■'9(WP<riwwMrii^^  isuouo  Mmoo’  29000^ 


fSTATCMCtTpi 
.  ffUMSER 

[dniT 


111  1 1  > 
Ml  I  I  ' 


|ii.| 

Odd  ■ '  ■  0  00 

I  I  ■  nil  11  n  iin  n n  ii  nu  i 

III  I  1  11  I  I  I  '  . . 


FORTRAN  STATEMENT 


dHT 

:i  11  I*  ; 


0  9  0  ' ' '  d  0  d  ’  '  III  '  '  D I . Odd  '  Td  ' 

n}Tifn3iit82ssiX3i}7a3«4i0  4i«e«<i«ossius3js»st»usitcc'Jk:i'<o»t:tt.a 

'III  III''  •  •  •  • '  1 1 1  n  1 1 1 1 ' ' '  •  ■  ■ 


:  f  KTir  I'.aTum 


oeouiTmToi 

tv  WTaN  nn  JTTina^ 


Frequencies  =  3.0000,  5.  0000  ....  30.0000 

Note:  This  frequency  card  should  contain  only  frequencies  of  the  high 
frequency  band  (3tp30  Me/  s)  unless  fewer  than  ten  frequencies 
are  desired.  In  this  ease,  enter  990000  in  the  unused  fields 
at  end  of  complement. 

Figure  IZ 


[Br  npsri 
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Example  Method  Card 


lwlf 

Foa  2 
VcOBMrtNf  * 

STATtMCNl ! 
NUMSER  ! 

1  -  . - . . . . .  -  -  ■  ---  - 

IDENTIFICATION 

■  fj^TRAN  STATE  ME  N  T 

unnrin 

lit  i  <  s 

r  e  0  0  3  0  fi  o'8  a  g  a  0  fl  11  mi  0  mi  i  ii  mil  iii  nos  (rim  0  DTiT'roTDd  d'o'  0  d  0  Ff  e  pi  u  0 

1  1  1  Hit  tt  DMIS  lltni  tIMlIZI  2aHflSnail»illl1ilia«])9l]|«4l4}0  4l«4|4l4l4t8bSI  S2U»«9SS4S7MS|{0tl  62I}U49H  47  64  M7C7I  U 

OOOOMtO 

TJ  14  19  14  11  71  71  M 

Example  Circuit  Card 


1 

ro« 

STATEbsEnT 

NUMSER 

wrr 

1 13  a  4  9 

1 

1 

1 

1 

fortran  STATEMENI 

IDCKTIFICATION 

InUT  irTTiTTOTr  oQeeiisii>ioas8fioooa  Tiriry  OoaaoootoftoaooooooTooT 

7  1  1  Ull  l>i3MI3l4tll|t«JiII  22  2)a4aJi9  2tiiatlAaia41l»a}ai]l49  4l041«4  484(47  4i4IN9l919>M»S4STMSIMIII2aM«lil7illl7en7: 

ooifinnriig 

n  74  79  71 77  74  79  tt 

Transmitter  latitude  =  46.  10  North 
Transmitter  longitude  =  64.80  West 
Receiver  latitude  =  32.  30  South 
Receiver  longitude  =  64.80  East 
Transmitter  power  =  1  kilowatt 
Transmitter  antenna  height  =  25  meters 
Transmitter  antenna  leg  length  =  0 
Transmitter  antenna  tilt  angle  =  0 
Transmitter  antenna  type  =  •£  (vertical) 

Receiver  antenna  height  =  30  meters 
Receiver  antenna  leg  length  =  180  meters 
Receiver  antenna  tilt  angle  =  70  degrees 
Receiver  antenna  type  =  -I  (rhombic) 

Required  S/  N  =  42  decibels 

Bearing  of  transmitting  rhombic  =  0  (on  great  circle  or  not  rhombic  type) 
Bearing  of  receiving  rhombic  =  0  (on  great  circle  or  not  rhombic  type) 
Abbreviation  for  receiving  location  =  BD 

Figure  13 


Tb-  37.053 
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EXAM'Fl^E  DATA  DECK  FOR  RUN  UNDER  TWO  METHODS  FROM  CARD 


:ll  _  6  3" 

p9993j9p9999999999Si9 
I  153^-112058  2295-12020 
i|  38971  I  7468  3232  6473 


to  .32250 '76-1  30180  67-t  -3  48 
1  15  ■,  -2  66  96  67-1  -2  61 


50000  70000  90000  ItOOOO  130000  iSOOOO  2O0O00  1500DO  300000 


-1 

12  53 
It  60 

999999999999 


I  153^12050  3401  -7165  20  15  -3  25  96  64-1  -3  5!  Pl^ 

I  31^  (-7200  i520--t2058  40  85180  67-1  36108  66-1  -4  49  102  25 

I  340l|  |-7l65  1530  -3900  tO  81201  70-1-12  -2  -3  47  fiflj 

300pp  50000  70000  90000  100000  120000  150000  170000  180000  250000  | 

_ _ j 

FORTRAN  3TATE1/ENT 

•IllMC  •«  ^  _ 

Wins,  jj  dFinnruirswsTiiimTiriciiiicjsiTJiirs c  jTfWToiii 05 1 1  r? ; 1 1 1 1:  i I’ci ; .  1 1 : 

III  1  «  1  ij:  *  *  on  » « »•  M  •>  s :  .i.’ijtna  1  *71 K  a  M  '  »  J  a  K  4^ 4«  ‘  V  •  '.v  ;«  a  •  i 

lu  ! ,  1 1  u  m  u  n  1 11 ;  I  III  n  I  n  J 1 11 1  M  1 1 1  n  1 1 1 1 J  J  n  ' :  n  1 1  i :  I !  1  M  1  j  I  •, !  M  1 


i!2  :  2  2  2i2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  ’  2  ?  2  2  2  2  2  2  2  :  2  2  2  2  2  2  2  2  ::  2  r  ?  2  ?  2  :•  1 

I 

3i2  3  3  3  313  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  J  3  J  3  3  3  3  3  3  3  3  3  5  3  3  3  3  3  3  5  3  3  5  3  2  S  2  2  3  2  3  3  3  3  J  2  •  3  i  2  S 

j'i  4  <  4  4  4  i4  44444l444444444444444444l444l4l4l4444444(<44i44i4444M</< 


:?  2  2  2,;i:  2  2212 


4  4  4  4,«  44  <  W  4  < 


S'G  6  8  41$  8  8  6  8  S  8  4  6  6  6  6  8  6  4  I  (  6  4  I  S  I  C  4  M  6  6  44  8 1 8  n  4  e  E  B  8  e  e  S  8  S  4  ^  U  E  {  r  :  E  E  U  f  (  k  [  2 :  U  [  H  t  :  r  £  t  i  1 
’ll*  '  i 

!  I  '  107 1  it!?  7  7  7  T  1  1  7  7  7  7  7  1  7  T  7  n  7  n  7  7  7  7  7  7  j  7  7  7  7  7  7  7  T  7  7  7  7  7  7  7  7  :  7  7  7  i  7  7  7  7  7  T  T  :  1 

•!)  19  ST. 1 1 1 1  9  •  1 1 1 1  •  •  •  i  1 1 M 1 1 1 1 1 1  •  1 1  •  1 1 1 1 1 1  •  9!  i  1 1 1  n  i  7 1  e  C  H  U  E  C  n  [  1 1  t  1  t  i  1  7  i  9  i  7  I  1  1  1  r 


Figure  15 


EXAMPLE  DATA  DECK  FOR  ONE  METHOD 
USING  CIRCUITS  FROM  MAGNETIC  TAPE 


1  iaitunMmiutiiiaiiiini<aaitiiaa)iiiiiiiaa»aB«ii«aitii«»aii«uuu>iuiis>iiaiaiiuu«iiiiiiiiun;i  dii»nKi"i»al 


1 1  n  1 1 11 1 1  n  1 1 1 1 1 1 1 1 1 1 1 1 1  n  1 1 1 1 1 11 1  n  I  n 

11112221222222222222222222222222:222222 

}2l2}22122123i3323223332332222223322233 
444444444444444444444444444444444444444 
53S$S555SSS355555S9iS5S3S33Si35555SSSSS 
till  4  4  6  441  1  6  4  11  4  4  (  6  41  1  1  t(lll(4(((S4SaeB 

I 

2  1  2  7  7  7  2  7  2  7  2  7  7  2  7  2  7  r  2  7  7  n  7  7  7  7  ?  7  7  2  7  2  7  7  3  7  )  7 

llll•l(l(ii((l((4((((((((((lil(((6(((((1 


,111111111111111 

11IM11IM 

222222222222222 

2222222222 

22212222 

3.33  3  333  3  3  3  3  33  3  3 

3333333333 

444444444444444 

4444444444 

44444444 

3S95SS339S54S9S 

9443333331 

(((((44(((((((( 

(((((((4(1 

777727777777777 

7777727772 

22222272 

((ll(•(((((((4( 

((1(311111 

Figure  16 
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EXAMPLE  HF  PREDICTION  DECK  TO  EXECUTE  UNDER  IBM- 7090  FORTRAN  MONITOR 


FORTRAN  STATEMENT 
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Vll.  PROGRAM  OUTPUT 

Tbe  output  is  showrt  in  Figures  i  througk  10.  Tbie  Output  is 
designed  for  9  by  11  standard  Gomputer  paperi 

A  description  of  variables  appearing  in  tbe  print-rOutS  reading 
left  to  rigbit: 

1.  Sednential  circuit  number  appearing  in  tbe  eidtreme  upper 
left  is  generated  within  tbe  macbine  and  is  designed  to  be  useful  for 
ready  reference  to  a  master  list.  Sample:  ll8 

li  Month  of  year  for  wMcb  predictions  are  calculated. 

Sample :  MARCH 

3.  Sunspot  number  used  in  predictions*  Sampe:  SSN  =  10* 

4.  Customer  identification  number  which  is  generated  from  a 
two«letter  code  for  receiver  location,  bearing  (receiver  to  transmitter)t 
and  great  circle  distance  of  path.  Sample:  KD  12.  012 

5.  Transmitter  geographic  coordinates  in  hundredths  of 
degrees.  Sample:  13.4SN  -  144.  80E 

6.  Receiver  geographic  coordinates  in  hundredths  of  degrees* 
Sample:  26.  72N  -  127.  80E 

7.  Great  circle  bearing  of  transmitter  to  receiver  and  bearing 
of  receiver  to  transmitter.  Sample:  312.4  •  126.5 

8.  Type  of  antenna  used  at  transmitting  terminal.  Sample: 
Rhombic  12H  39L  58  Deg. 

H  -  height  in  meters 
L  leg  length  in  meters 
Deg  =  tilt  angle  in  degrees 

9*  Level  of  man»made  noise  at  receiver  site.  Sample: 

NOISE  =  3 

10,  Type  of  antenna  used  at  receiving  terminal.  Sample: 
VERTICAL  *  14H  OL  ODEG. 


-  45  - 

The  al4H  indicates  the  grounded  vertical  is  1/4  (\).  The  sanie 
analogy  would  apply  for  -12H{X./2)i  «iH  (\),  etc.  A  print »out  without 
the  minus  indicates  the  vertical  height  in  meters  (e.g. ,  i3H  indicates 
the  grounded  vertical  is  13  meters  high). 

11.  Great  circle  distance  of  path  in  nautical  miles.  Sample: 

N.  Miles  1243.2 

12.  The  effective  power  delivered  to  the  transmitting  antenna. 
Sample:  PWR  s  20.00  kw 

13.  The  required  signal -to -noise  in  a  one -cycle  band  (db  >  1  watt) 
at  receiving  antenna  terminals.  Sample:  REQ.  S/ N  =  45  db 

14.  Greenwich  mean  time.  Sample:  GMT 

15*  Maximum  Usable  Frequency.  Sample:  MUF 

16*  Optimum  Traffic  Frequency.  Sample:  FOT 

17.  Truncated  numbers  representing  the  frequency  complement 
card  (Figure  14).  Sample  4,  6,  8  . . . .  23 

18.  Probable  mode  associated  with  least  loss.  Sample:  3Fi  4E 
and  4X.  4X  indicates  propagation  by  E  and  F^  region  using  four  hops. 

19*  The  theoretical  angle  of  arrival  associated  with  the  most 
probable  mode  of  propagation.  Sample:  20 

20.  Reliability  defined  as  the  per  cent  of  days  with  the  month 
the  available  signal-to-noise  will  be  adequate.  Sample:  Reliability 

21.  Signal-to-noise  defined  as  the  available  monthly  median. 
Sample:  S/N  .  .  .  .  db 

22.  System  loss  defined  as  the  ratio  of  the  power  delivered  to 
the  transmitting  anteima  relative  to  that  available  at  the  receiving 
antenna  terminals.  Transmission  loss  may  be  obtained  by  using 
isotropic  antennas.  Sample:  Loss  ,  .  ,  db 

23.  The  Lowest  Useful  High  Frequency  defined  as  the  lowest 
frequency  that  will  give  90%  reliability  at  a  given  hour.  Sample;  LUF 
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24.  A  zero  appearing  in  the  LUF  column  denotes  a  LtJF  greater 
than  the  FOT. 

ZBi  A  a3  in  the  LUF  column  denotes  a  LtJF  below  3  Me/  s. 

26.  A  calculation  in  any  method  showing  all  zeros  denotes  a 
frequency  for  which  the  calculation  was  not  completed  due  to  mode 
limitations  or  the  primary  calculation  for  the  desired  mode  exceeded 
the  absolute  value  999^ 
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Vni.  COMPLETE  FORTRAN  LISTING  OF  COMPUTER  PROGRAM 
FOR  PREDICTING  MF  SYSTEM  PERFORMANCE 

C  MUfLUF.**«H<^F  PREDICTIONS  BY  D«L«LUCAS  AND  6*W*HAYD0N  ••••••CRPL..... 

C  FORMATS  AND  STORAGf  ASSlGNHtNtS 

6  format (10F7 .4)  0001 

1  F0RMAT(2(FS^2iF6.2) . I2t5>F3.0»I2*3F3.0*I2*I3fF3.0»2F4.0i4XtA2)  0002 

3  F0RMAT(10X*I5*10X*A6»  lOX»4HSSN»FS^0*10X*A2#F7.3/9Xt  UhTR  0003 

2  ANSM 1 1 1 ER • 1 3  X  «  SHR  ECE I VfR  *  1 2X 1 8HBE AR I NGS  *  6X 1 7HN • M I L  E  S / 7X  « F6 • 2  *  A 1  *  0004 

32M  «  •F7.2*AltSXtF6*2*Al*2H  ^  •F7.2*Al*3Xt2F6»l*4XfF8.I )  0005 

4  format  (FIO.O)  0006 

5  FORMAT (41 2)  0007 

6t6  FORMAT! 12A6)  0008 

7  FORMAT (lX/6Xt  0009 

43HGMT  * 2X • 3hMUF  *  3x  * 3hFOT  * 3 ( 5X  »  3mGMT  *  2X  » 3HMUF  » 3X  »  3HF0T ) /  IH04( 4X  *  1 3  *  0010 

52F6.1)/lH04(4X*I3*2F6«i)/lH04<4X»i3»2F6.1)/lH04(4X» I3*2F6.1>/1m04(  0011 

64X«I3t2F6.1)/lH04(4X«I3t2F6*l)/lH  Al)  0012 

DIMENSION  AM0N(12)>  GLAT (5 ) >ABI ( 24 ) • ABI Y i 5 »24 I •CLCK t 24) *07(5 ) *  0013 

1F2S(5*24)iF2H(5*24)*GML(S*24)*6MM(5*24)*MF(3)»ME(2)  0014 

2*ELO(10) »FLD(10) tF24(5 )*GMA(S)»ADJ(4) •ABC( 12)  0015 

DIMENSION  RD(5) *SUN(12)*A(10*7»l4) »RaSSN( 12 >  >EKF( S»24)  OOU 

DIMENSION  0(20*60t4)*IL(4)»JL(4)>KU(4)»LX(4)»JAL(4)  •CKC(24)  0017 

1 fG(60 ) tAB(60 ) »S( 20*24 ) tC( 20f24 )»GAMMA( 5*24 )»BA( SO)  0018 

DIMENSION  CLAT(5)*CLONG(5>*EMMFY(3t24)*FMUFY(&»24)»  UFY ( 5 ) »FOTY ( 5  )  0019 

5*16(24)  0020 

PIMENSION  UF(24)*F0T(24)*P(29*16*5)*ABP(2»6)  0021 

dimension  AHA(2) *ARA(2)*AVA(2)  0022 

DIMENSION  X1XX(12)  0023 

DIMENSION  FREL(12)  0024 

COMMON  ABP*P  0025 

COMMON  A  *9  *  G  *  AB  *  S  *  C  *  G AMMA  *BA • EMUF Y . FMUFY  *  I G  tUF  *  FOT • SUN  *  AVA  *  ARA  *  QU A  0026 

1T*ABI*ABIY*CLCK.*GT*F2S»F2H*GML*GMM*MF*ME*ELD*FLD*F24*GMA«ADU*ABC*  0027 

2EMF*CKC*CUAT*CL0NG*AK*BX*CX*pX*Pl2*EX*EEX*GLT*GLG*AAA*ALA*ASA*  0028 

3AFC*ANC*AWC*AEC»ASC*A2Z*MaN*RSN*XM>PWR*|ANT>IANR*Y2*X2*MOUSE*SSN*  0029 

4N0CIR*K4*  AX*XLONG*AU>YtAT*XLAT*YLONG*YL.*BTRY*GCDXM*GCDNM*Xl*K2*K  0030 

53*K5*GCD*rRSN*MON*ID*MOUSE*AY •METHOD* AHA*BM0NS*MAP  0031 

COMMON  FREL  0032 

COMMON  XNH*XNL*XND*RNH*RNL*RND*HA*JIG*HAR*XTR»XETA*WjETA»IHR*MIT  0033 

REWIND  2  0034 

REWIND  3  0035 

909  format (IMl)  0036 

DO  101  I«l*24  0037 

101  I6(I)«I  PP38 
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C 


B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 

6 

B 

B 

B 

B 

B 

B 

B 

B 

B 

B 


CONSTANTS  TO  BE  USED  PREOUENTLY 

AK«l*7A5329E-2 

BK«5»729577El 

CK>i.lll36E2 

DM  .0062 137 

PI2-1.37079 

Exa. 539956 

EEK«.6214 

6l.T«78i*AK 

GL6>70«»AX 

ANON ( 1 1 -412145666060 

AMON ( 2 ) -262522606060 

ANON ( 3 ) -442151606060 

AMON (4  >  >214751606060 

AMON ( 5 ) -442 1 70606060 

AMON ( 6 ) -416445606060 

AMON ( 7 1 -416443606060 

AMON (8 ) -216427606060 

AMON ( 9 ) -622547606060 

AMON ( 10 ) -462363606060 

AMON (11) -454665606060 

AMON ( 12 ) >242523606060 

AAA-606060606060 

APC-266060606060 

AMA( 1 ) >606060603040 

AHA(2)«243147464325 

AVA( 1 ) -696060606525 

AVA( 2) -516331232 143 

ARA  < 1 ) >606060606051 

ARA( 2) >304644223123 

AZ2>006060606060 

ANC -456060606060 

AWC-666060606060 

AEC-256060606960 

ASC-626060606060 

PNINE-llllllllllll 

SUN(1)>«21. 

SUN(2)-»13. 

SUN(3)-«03. 

SUN(4)-09« 


0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

005B 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0060 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 


« 
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SUN(5)»18.  6678 

SUN(6)>23.  6679 

SUN(7)«21.  6680 

SUN (8) ■14.  6681 

SUN<9)>64.  6682 

SUN<10»«=68,  6683 

SUN(il)ii^l4.  6684 

SUN(12)>^23.  6683 

ADJ(il*9.  6686 

ADJ(2)«174  6087 

AOJ(3)B26*  6688 

ADJ(4)a28.  6689 

C  READ  CIRCUIT  DATA  AND  SEARCM  DATA  TAPE 

WRITS  OUTPUT  Tape  6 *909  0096 

464  READ  INPUT  TAPE  3*3*  6691 

1  METHOD*METH*IhR*MAP  0092 

READ  INPUT  TAPE  3t6,(FRELM)*I*l*l6l  0693 

Mlful  0094 

IF  (METHOD!  463*12*463  6695 

12  Rewind  2  0696 

WRITE  OUTPUT  TAPE  6*909  6097 

Rewind  3  6098 

CAUL  Exit  6099 

463  Ll>2  6106 

60  TO  l30*30*30*3d*3l*30)«METMO0  0101 

30  Ll«l  0102 

3|  N0CIR>  0  0163 

MONTMaQ  0104 

IFIMETH!  7770*7771*7770  0103 

7770  READ  INPUT  TAPE  3*676*XlXX  6166 

WRITE  OUTPUT  Tape  3*676»X1XX  6107 

IFIFNINE’XlXXdl )  7770*7771*7770  0108 

7771  REWIND  3  0109 

83  read  input  tape  3*1*  0116 

1  XiX*YlY*X2X*Y2Y*MONS»SUN5*TWR*XNHO*XNLO*XNDO*ITANT*RNMO*RN  Olll 

1 LO  *  RNDO  *IRANT*MANN*RRSN*XETA*WETA*HAP  0112 

16  XlvXlX  0113 

Y1»Y1Y  0114 

X2*X2X  0113 

Y2«Y2Y  0116 
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IANf«ItANT  0117 

lANR-IRANt  0118 

MAN -MANN  0119 

RSN-RRSN  0120 

iRsN«RSN  oiai 

PWR-TNR  0122 

XNN«XNHO  0123 

XNL>XNLO  0124 

XNOiiXNDd  0125 

RNHiiRNHO  0126 

RNL-RNLO  0127 

RNOaRNOO  0128 

HAaHAO  0129 

iF(MdNS)  464ait82  0130 

02  LOCKS  MONS^HdNtM‘1  0131 

SSN-SUNS  0132 

MONTHsHONS  0133 

MOUSEsMONS  0134 

KZ>1  0135 

READ  INPUf  tAPE3*l*XlX«YlY*X2X*Y2Y>MONS*SUNS»twR*XNHO»XNLO»XNpO»  0136 

1 ItANT*RNMOiRNLO*RNOQt iRANt •HANN*RR5N*XEtA»WETArHA0  0137 

GO  TO  10  0138 

11  GO  to  ( 15*77 )*XZ  0139 

15  NOCIRSO  0140 

WRITE  OUTPUT  TAPE  6*909  0141 

62  JiG>0  0142 

MIT«1  0143 

KZs2  0144 

IF(LOCK)  77*731*33  0165 

33  00730  IU>l*LOCK  0146 

739  READ  TAPE  2  0147 

731  read  tape  2*IL*JL*KL*LK*JAL*Q*A*P*ABP  0148 

77  READ  INPUT  tAPE3*l»XlX*YlY*X2X.Y2Y*MONG*SUNS*TWR*XNHO*XNLP*XNpO*  0149 

1 1 T ANT  * RNMO  * RNLO  *  RNpp  » 1 R AN T  *MANN  * RRSN  *  XE T  A  *W | T  A  t  MAO  0150 

IF(91.-XIX)  78*78*16  0151 

78  REWINP  3  0152 

READ  INPUT  TAPE  5*1*  0153 

1  XlX«YlY*X2X*Y2Y*MONS*SUNS*TWR*XNMQ*XNLQ*XNpO* ITANT*RNH0*RN  0154 

1LO*RNPO*IRANT*MANN*RRSN*XETA*WETA*HAO  0155 

IF(MONS)  83*16*83  0156 
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83  IFIMONS-MONTWI  84 •16.16  0157 

84  REWIND  2  0158 

MONTHS Q  0159 

16  SSPSSUN (MOUSE )*AK  0160 

BMONSSAMONI MOUSE )  0161 

NOGlRsNOGlR+l  0162 

J1P*1  0163 

JEPsl  0164 

G  GREAf  CIRCLE  DISTANCE  AND  DISTANGE  TO  CONTROL  POINTS 

XXl#Xi»AX  0165 

XX2sX2*AK  0166 

YYlsYlSAK  0167 

YY2SY2SAK  0188 

GCOsSI  NF  t  XXl » *Sl  NF  ( XX2  l+COSF  ( XXI )  *G0SF  ( XX2  )  *G0SF  (  Y  Y2“-YY  1 )  0169 

GGDsAGOSF(GGD»  0170 

RD(1)-GCD/2«  0171 

RD(2)«18**AX  0172 

R0(3)«GG0-18.*AX  0173 

RD(4)s9**AX  0174 

RP(5)«6G0“9.*AK  0175 

GCOXM*GCD*eK*BX  0176 

6C0NMsGCDXM*EX  0177 

C  SEARtNGSS-”^  FORWARD  AND  BACKWARD 

OFSYY2-YY1  0178 

IF(ABSF(DF1-3*141592654)  26*22»22  0179 

22  IF(PF)  24*25i25  0180 

24  OFS6.283185308+DF  0181 

GO  TO  26  0182 

25  DFi-6«2831853G8*OF  0183 

26  y»f5*(Pl2’XXl*Pl2-'XX2*GCP»  0184 

eiN0«SINF(U)»SINF(U*PI2*XXl)  0185 

IF(eiND)  27.28.27  0186 

27  BTRVsll4.5816*ATANF(SQRTF(ABSF(SlNF(y-PI2+XX2)*SINF(U-GeD) »/eiNO) »  0187 

IF(OF)  7029.21.39  0188 

21  IF(XX1^XX2)  28»39.39  Ql89 

28  BTRY-IBO,  0100 

GO  TO  39  0191 

7029  BTRY«360.’BTRY  0192 

39  RF8YY1»YY2  0193 

IF(ABSF(RF)*3, 141592654)  126.122.122  019^ 
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122  ifCRF)  124tl25*125  019S 

124  RFa6»283185306‘^RF  0196 

GO  TO  126  0197 

125  RF»*6* 283185 368+RF  0198 

126  Ua*5»(P12*XX2+P12-XXl+GCD»  0199 

C IND  a  S 1 N  F  t  U ) *S I N  F ( y «  P 1 2 + X  X  2 »  0  200 

IF(CINO)  57.778.57  0201 

57  XTRaii4,5816*ATANF(S0RTF(ABSFiSlNF(U-Pi2‘^XXl)*SiNF«U-GG0) »/CiNO»)  0202 

IF«RF»  8029*51.79  0203 

51  1F(XX2»‘XXI)778*79»79  0204 

778  XTRalSO.  0205 

GO  TO  79  0206 

8029  XTR  a360.-xTR  0207 

79  iFtAOOO.-GCOKMl  40. 41. 41  0208 

40  Xlai  02O9 

K2a2  0210 

X3a3  0211 

X4a4  02 i2 

<5*5  0213 

GO  TO  43  0214 

41  Rial  0215 

<2*1  0216 

<5*1  0217 

<4*1  0218 

<5*1  0219 

C  CONTROL  POINT  LOCAL  TIME  .  G|OGRAPHIC  LATITUDE  AND  LONGiTUDE 

43  DO  61  ITaL1.24.Ll  O220 

GMT  a  IT  0221 

G0LOaGMT-Y2/l5.  0222 

IF(24.«^COLO)  94.95.95  0223 

94  COLOaCOLO«24.  0224 

GO  TO  96  0225 

95  IF (COLO)  98.98.96  0226 

98  C0LO«C0L0.>24f  0227 

96  C<C(IT)a0OLO  0228 

8SLa(15.«GMT-180.)  0229 

DO  61  I|aKl.K5  0230 

PPaRDdl)  0231 

CENLATaCOSF  ( PP )  •COSF  ( P 1 2*XX2 )  I NF  ( PP I  0232 

l*SI NF ( P I 2-XX2 ) tCOSF ( BTRY*AR )  0233 
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CfNLAT«(PI2-AC0Sf (CENLAT) l*BK  0234 

GENLG*  i  KGSf (PP)*COSF»Pi2-XX2)*GGSF(Pl2-GiNtAT*AK) >/<SiNF(PI2  0235 

i-XX2)*SiNF(Pl2^GENIkAT*AK» »  )  0236 

lF(ABSF(GiNLG)»UI  3001  *3001  ♦3000  0237 

3000  GENLe*!#  0230 

3001  GENLG«AGOSF(GENLG)  0239 

GENLG*  <  YY2^S IGNF ( GENUQ.DF ) )*BK  0240 

1F<180.»ABSF(CENLGJ )  7it69*69  0241 

71  IFIGENLG)  73»72*72  0242 

73  CENL6»360,-*CENLG  0243 

GO  to  69  0244 

72  CENLGtGENLG-360*  0245 

G  GEOMAGNEtiG  LOCAflON  OF  GONtROL  POINtS 

69  6At»  AGOSF  (  Sf  NF  (GLt  t  *5  INF  I  GENtAT»<kX  )  ♦COSF  ( GCf )  *GOSF  ( GENLAt»AK  )  *  0246 

IGOSF  <  gENL6*AK^GLG ) J  0247 

GUAT< ii )«PI2*GAf  0246 

GLATdnsGENLAT  0249 

70  GLON6( li»*GENUG  0260 

C  SUNS  ZENITH  ANGIlE 

Z«(SSU*G|NLG>*AK  0251 

56  GYG6N«SINF  <  GENLAt*AX » *Sl NF  <  SSP ) ♦COSF  <  GENt AT*AX ) *GOSF  C  SSP » *CGSF I Z »  0252 

GyGEN«AGOSF<GYGEN1  0253 

GYGtN*AB$F(GYGENl  0254 

C  ABSORPTION  INDEX  I  AT  CONTROL  POINTS 

32  1F( l02*tAK*GYGEN)  50t58f59  0255 

58  ABIY<II*IT)>0.  0256 

GO  TO  886  0257 

59  ABIYnitlTl.(l,*,0037#SSN)*GOSF<,881*GYCEN)«*l,3  0258 

886  GO  TO  (165tl66)tK2  0259 

165  ABI(tT)>ABIY(IItIT]  0260 

GO  to  61  0261 

166  ABK  IT)-(ABIY(lflT)*^ABIY(2tIT)4-ABlYI3flT)  )/3.  0262 

61  CONTINUE  0263 

DO  2112  II-KltK5  0264 

DO  2112  lT>Llf24fLl  0265 

AB|C>ABIY( II flT)  0266 

C  EHAYER  PI  STANCE  FACTOR 

ARC>PK#GCPRM  0267 

1F(16*^ARC)  88*87f87  0268 

87  ELFG-I  <  ( ( I  («4.368460907E’'9*ARC-*-1«334494261E^7)*ARC’-5»976618436E«6I  0269 
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l»ARG+2.624808315E-4)*ARC-5*038476266€-3)*ARG-t»3,761385053E*2»*ARG-l  0270 

2*1332 007  56 E “ 2 ) « AR  G  * .2085  0271 

GO  TO  100  0272 

88  ELFG«1*02  0273 

G  E-2000  CONTROL  POINT  MUF 

100  |MG«( ( (-iO,66484998*A6iG+39.26l5lO56)*ABlC-52.41i0l754)*ABlG+37,67  0274 

1726072 )*ABlG+3, 345996232  0275 

EMF<  niiTl»EMG  0276 

2112  iMUFY« I I*iT)sEMG*ELFG  0277 

G  GENERATE  CONTROL  POINT  F0F2  ANO  M-30G0  FACTORS 

00  1112  iI-K2«K3  0278 

504  GENLG«GlONG( 11)  0279 

GENLAT-GLAT( 11 )  0280 

IF(GENLG)  4l2t4l4t414  0281 

414  GLGs360.^GENLG  0282 

GO  TO  413  0283 

412  GLG«AB$F(GENLG)  0284 

413  B0Y>SINF(GLG*AX)  0285 

C0G«SINF(2.*GLG*AK)  0286 

OOG«GOSF(2«*GL6*AK)  0287 

G06«GOSF(CENLAT*AX)  0288 

606«$lNF(GgNLAT*AK)  0289 

H0G>C0SF(CLG*AR)  0290 

DO  500  I0«lt4  0291 

415  I*IL(I0)  0292  ‘ 

J«JL ( I 0 )  0293 

KvRLdO)  0294 

LsLK(IO)  0295 

JAipJALnO)  0296 

DO  908  KAalrK  0297 

KPaXA"!  0298 

408  G(KA)-BPG9»KP  0299 

LA»0  0300 

tP-L-l  0301 

XRvK*’!  0302 

DO  409  KA-KXtLO*2  0303 

G(XA)»B009*LA«QPB«HOC  0309 

6 ( KA41 ) -BOG«*LA*G0B»B0Y  0305 

909  LA«LA+1  0306 

LB«0  0307 
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0308 

0309 

DO  410  KA>lLL*iMi2 

0310 

G  <  KA ) >B0G4*4B«G0B*«2»00G 

0311 

GIKAi-l )  -b0g*«4B*GOB**2*COG 

0312 

410 

LB-LB-^I 

0313 

DO  411  JBai*J 

0314 

ABIJB)aO. 

0315 

iS-2«JB»i 

0316 

DO  411  KA>l«i 

0317 

411 

AB <  JB ) a AB 1 JB ) *Qi  1 5 *XA * lO t  *G 1 KA ) 

0318 

DO  407  jBa2*J 

0319 

BA(JB)aO« 

0320 

IS»24JB*2 

0321 

DO  407  KA-UI 

0322 

407 

BA( JB)aBA<JB)aQ< IS«KA*I0I*6(KAI 

0323 

DO  500  lt>4lt24.Li 

0324 

503 

QMtalf 

0325 

CLOCKaGMt-CEN46/l5, 

0326 

IP(24*^C40CK)  64*65«45 

0327 

64 

Gt.OCKaCLOCK^24. 

0328 

GO  TO  66 

0329 

65 

IF (CLOCK)  68«6e*66 

0330 

68 

CL0CKaCL0CK-f24« 

0331 

66 

CLCK(IT)«CL0CK 

0332 

TIMEa(19,ttCLOCK»180.)aAK 

0333 

DO  998  JB>2*J 

0334 

FBaJB^l 

0335 

S(JB»IT)»SINF(FB*TIME) 

0336 

998 

CIJB>|T)-C05F(FB*TIME) 

0337 

GAMMA! I0*IT)-AB(1> 

0339 

DO  500  JB-2tJ 

0339 

500 

GAMMA! lOt IT) -GAMMA! lO* IT )4>AB! JB)«C! ^B* 1 T laOA! JB)aS! JB t I T ) 

0340 

c 

F2’^REGIQN  GYRO^FREQUENCY 

152 

CALL  POLY ! 1 f 7 1 7 # CENLQ/ 180. »CCNL AT /90, » Y ) 

0341 

GYR0aY/2f 

0342 

GYdDaY 

0343 

C 

F2^LAYER  CONTROL  POINT  F2o4000  MUF 

DO  1112  ITaLl«24,Ll 

0344 

F2LS-GAMMA ! 1 1 1 T I aGAMNA ! 3 1 1 T ) • 1 . 1 

0345 

5  6  * 

F2StIi#IT)»F2US  6346 

f2HS«GAMMA ( 2 • I T ) «6AMHA ( 4  *  I f ) *1 ^ 1  6347 

F2M< li *iT)>F2HS  6348 

G  F2«LAYER  OlSTANGi  FACTOR 

lFt24*-ARC)  157.157*158  6349 

157  FLFGsl.  6356 

60  to  159  6351 

158  F6FCa( ( ( ( < (*6.7l2654756E“9*ARG»4*4915l44lt-7)*ARG-9.98583ll64l«6)*  6352 

iARG+6. 8652598171*5 )*ARG+9 .2624373328*5 )*aRG‘»2. 2646343411*3  >»ARG+4.  6353 

26992431618*3 )»aRG  6354 

159  GAHMAdiif  )a6AMMA(l.it)4-6YR0  6355 

GAHMA ( 2  *  1 1 ) aQAHMA ( 2 • I T ) ♦GYRO  6356 

FMLSsGAMMA(l»IT)+FLFG*(F2US*GAMMA«l»it))  6357 

GMU ( t i * i f ) >GAMHA ( 1 « 1 1 )  6  358 

FMHS»6AMMA  <  2 . i t ) ♦FlFG* ( F2HS*GAMMA ( 2  *  1 1  n  6359 

GMH( II •iT)aGAMMA(2.lt|  6366 

C  F2*tAY8R  CONfROL  FOINt  MUF 

FMUFYi 1 1 tiT)a(FMUS»< i86*-SSN)*FMHS*16SN*10.) )/176*  6361 

1112  CONTINUE  6362 

G  G I Reg it  MUF  ANO  FOT 

00  1612  IT>L1*24.^1  6363 

6MT«IT  6364 

60  1060  II#K2*K3  6365 

IF(K2^1I  662.663*662  6356 

663  |R«1  6367 

GO  TO  664  0368 

662  I Rail ♦2  6369 

604  |F(EMUFY(IR.|T)*FMUFY< II.IT)»  1007.1007.1008  6370 

1007  UFYIinaFMUFYdl.IT)  6371 

GO  TO  1009  0372 

1008  UFY( II )aEMUFY( IR.IT)  0373 

1069  F0TFYa,85*FMUFY(II.IT)  6374 

IFdiMUFYUR.ITI-FOTFY))  lOlO.lOll.lQll  0375 

1610  FOTYI  inaFQTFY  0376 

GO  TO  1000  6377 

lOll  FOTYIIIIaEMUFY(IR.IT)  0378 

1060  CONTINUE  0379 

lIaK2  0380 

IFtl^in  606.605.606  0381 

605  UFtlTlaUFYIl)  0382 


704  AY»A$C  0409 

705  YL«AWC  Q410 

IF(Y2)  7Q6*2221t2221  0411 

706  YL*AEC  0412 

2221  GO  to  (2223t2224*2224t2224>2224*22?4)*MEtH0P  0413 

2223  IF(MAP)  3223t3223t3224  0414 

3224  CALL  CURVY (FOT (OF)  0^15 

60  TO  IQ  0416 

3223  WRITE  OUTPUT  TaPE  6»3tNOClRtBMONS»SSN*HA>HARtXLAT i AX.XLONGtAL*  0417 

1YLAT»AY*YL0N6«YL«XTR»BTRY*6C0NM  0618 

WRITE  OUTPUT  TAPE  6t7.  0619 

1  (16(1)  »UF( ! )  tFOT(I )  tI6(  I’^6}  >UF ( IH >  »F0T  ( I-t<6)  1 16(  I<fl2)  *  0620 

2  UF(  I-i<12)*F0T(  I4'i2)tI6(  HIE)  tUF(  14^18)  »fOT(  14.18)  tI>l>6)*AAA  0621 

IF(3«MIT)  969»969*970  0622 
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969  WRITE  OUTPUT  TAPE  6.909  0423 

MIT«1  0424 

GO  TO  10  0425 

970  MIT«M1T+1  0426 

GO  TO  10  0427 

2224  gall  LUFFT  0426 

GO  TO  10  0429 

END  0430 
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SUBROUTINE  LUFFY  0001 

C  FORMATS  AND  STORAGE  ASSIGNMENTS 

01  MENS  I  ON  FDEKdOliEOEKdO)  0002 

dimension  AM0Nd2)t  6LAT(5)  •ABI  <  24)  •  ABIYC  5 124)  tCLOKI  24)  »GY(  5  )  t  0003 

1 F  2  S ( 9 1 2  A ) *  F  2H ( S 1 24 ) «  GML (9*24) »  GMH ( S • 2  4 ) «  MF ( 3 ) *  M  E ( 2 ) *E  L  L  d  0 ) *  F  UU (1 0  0004 

2 )  tELD dO )  tEDEL ( 10 )  »FUD  ( 10 )  *FiDEL  (10 )  * F 24 ( 3 )  *GMA  ( 9 )  * FUF (10 )  *FSK YdO )  0009 

3*ESKY(  10)  *GE)iO)  *GFdO).FsE(  lO)  rFSFdO)  *NOS(  11 )  .FXMT(  10)  *EXMt(10)  •  0006 

4FRCR  d  0 ) t  E  RC  R (10)*  ADU ( 4 ) *  ABC ( 1 2 ) • MOD  E  d  2 )  0007 

DIMENSION  Q(20*60*4)*IU(4)*JI.(4)*KL(4)*4XI4)*JA4(4)  *CKC(24)  0008 

1 *G(60 ) *AB( 60 ) *S( 20*24 ) *C ( 20*24 ) *GAMMa (9*24) (BA ( 60 )  0009 

DIMENSION  RD(9)*SUNd2)*AdO*7*14)  *RaSSN(  12  )  iEMFI  9*24  )  0010 

DIMENSION  GLAT(9)*Ct.0N6(5)*EMUFY(9*24)*FMOFY(9*24)*  UFY ( 9 ) *FOTY( 9 j  0011 

9*10(24)  0012 

D I MENS  I  ON  UF ( 24 ) * FOT ( 24 ) *  FF ( 24 ) * P ( 2 9  *  1 6  *  6 ) *  ABP ( 2  •  6 )  OO 1 3 

DlMENSiON  AHA(2)*ARA(2)*AVA(2)»RANT(2)*tANT(2)  0014 

DIMENSI0NNXL0S(  1 1 )  *NAN6(^E  ( 11 )  0019 

DIMENSION  TOP (4)  0016 

DIMENSION  IFRdO)  0017 

DIMENSION  FR|l.d2)  0018 

COMMON  ABP*P  0019 

COMMON  A*Q*G*A6*S*C*GAMMA*BA*EMUFY»FMUFY*lG*UF*F0T»SUN*AyA*ARA*GLA  0020 
1  f  * ABI  *  AB  I  Y  *ClCK  *  GY  t  F2S  *  FEh  *GM(^  *  GMM  .  MF  *  M| .  E  LO  *F  I.D*  F24  *  GMA  *  ADU  »  ABC »  002 1 

2EMF*CKC*CLAT*CL0NG*AF*BK*CK*DA*PI2*EX*EEK*GLT*GLG*AAA*ALA*ASA*  0022 

3AFC*ANC*AWC*AEC*ASC*AZZ*MAN*RSN*<W*PwR*lANT*)ANR*Y2*X2*MOUSt*SSN*  0023 

4N0CIR*K4*  AX*XLONG*ALtYLAT»XLAT*YLONG*YL*BTRY*GCDi5M.GCDNM*Xl*IC2*X  0024 

93  *R9  *  GCD  *  I  RSN  *M0N  *  1 0  *M01JSE  *  AY  *  ME  T  HDD  ,  AH  A  •  BMQNS  *  MA  P  0029 

COMMON  FREL  0026 

COMMON  xnh*xnl»xnd*rnh,rnl*rno*ha*jig*har*xtr*xeta*weta*ihR*mit  0027 

COMMON  i;AP*KUPtMA*TANTtlNM*|NL*INP*RANTtJNHtJNL*UND  0028 

2  FORMAT(2X*2A6*lX*I3ilHM*lX*I3*lHL*lXd3*3HPEG*12X*6HNOISEsI3*12X*4  0029 

1HANT-I3*2M0B)  0030 

1  F0RMAT(2X*2A6(lXt l3*lHH»lX*|3*lML*lXd3*3MDEG*4X*6MN0ISE*I3*2A6*lX  0031 

l*I3*lMM*lX*!3*lH(,*lXf  I3*3HOEG)  0032 

00  F0RMAT(6X*4MPWR"F6,2*2MXW*11X*21HQPERATING  FREQUENC|ES*9X»14HF|ELD  0033 

1  strength/ 9X*14HGMT  MUF  F0T*lX.l0I4t4H  FOT)  0034 

3  FORMAT(IOX*IS*10X*A6*  10X»4HSSN-F5,0*10X*A2*F7,3/9X»  UHTR  0039 

2ANSMITTER*13X.8HRECEIVER*12X*8HBEAR!N6S*6X.7hN,M!LES/7X  •F6»2*A1*  0036 

32H  *  *F7.2.Al*9XfF6,2»Al»2H  »F7,2 *Ai .3X*2F6. 1 .4X*F8,1 )  0037 

4  FORMAT(6X*4HANT»I3*2HPB»13X.6HNOISE«I3*7X*2A6*1X*I3*IHH»1X*I3*1HL,  0038 

11X*I3*3H0EG)  0039 


9  f0RMAT(4XtI3.2F6.1/20X.ll« I3tAl)t6H  MObE/ZOX* H j4t7H  AN6UE/4Xi  0040 

2l3t2F6*l.lXtliI4,l3M  RtUiABILItY)  0041 

12  FORMAT(lXtA2.F7i3,I5.2FlO»2l  OOA2 

13  FORMA  t ( 6X  r4  H AN  f  « I  3  «  2  MOB • 2  OX • 6HN0 1 S£  « i 3 • 2  0  X  »  4  H  A  N  T  >13*  2  HOB )  0043 

14  FORMAT ( 6X • 4HPWR« F6 . 2  *  2MKW • 1 1 X  * 2 ImOPIR AT  I NG  f RtOUE NC 1 1 S  » 1 OX . BhREO • B  0044 

1/N«13»2HDB/5X.14MGMT  HUF  FOT.1X»10I4*4H  FOT )  0045 

29  FORMAT(4X.l3t2F6.i/20x*il(I3»Al).6H  MODE/2GX*liI^*7M  ANGUt/^X*  0046 

2l3*2F6*l*lX*llI4t9H  5/N.,0&)  0047 

34  forma  T(6X*4HPWR>F6.2*2MKW»16X*1  3  H'R  i  L  I A  B  1 1. 1  T  1  E  S  » 1 3  X  »  8  M  RE  Q  t  S  /N  ■  1 3  *  2  M  0048 

10B/9X*14H6MT  MUF  FoT»1X*10T4»9H  FoT  MCS)  0049 

7  FORMAT  (6X*  0080 

43HGMT*2X*3MUUF*3X*3MFOt*3(5X.3>H6MT*2x»3HLUF*3x*3HFOT) /  iH04<4X*13*  0081 

52F6*l)/lH04(4X*i3t2F6^1)/lH04(4X*I3*2F6*l>/lH04(4X*l3*2F6*l)/lM04(  0052 

64Xtl3i2F6*l)/lH04(4X*13*2F6*l>/lH  Al)  0053 

8  FORMATC 31X*21MOPERaTiNG  FREQUENeiES/SX.lAHGMT  MUP  FOT»lX»iOl4*  0054 

14H  FOT)  0055 

19  FORMAt(4X*t3«2F6.i/20x*ll ( I3*A1)*6H  MO0E/20X*lli4*7M  ANGEE/4X*  0056 

213*2F6*1*1X*1114.10H  lOSS.*OB)  0057 

59  FORMAT(4X*i3*2F6.1/20x*ll(l3*Al)>6H  MOOI/20X* 11I4*7H  AN6tE/4X*  0058 

2l3*2F6«ltlX*llI4t5H  OBU)  0059 

24  FORMAt(6X*4MPwR«F6«2*2HKM»42X»8HREQ*S/N>t3*2MDB)  0060 

70  FORMAt<6X»3HGMT»2X»3HLUF»3X.3HFOT*3(  5X.3HGMT,2X*3MLUF*3X*3MFOT)/i  0061 

lH04(4X*i3*6X*F6«i)/lH04<4X*I3»2F6»l)/lN04(4X»I3»6X»P6*l)/lM04(4X»I  0062 

23'2F6*l)/iH04(4X*I3»6X*F6.1)/lH04(4X>t3*2F6*l)/lH  Al)  0063 

94  FORMAT(8X*i3*2F6.1tlX*lll4)  0064 

907  FORMAT< IHO/IHO/IHO)  0065 

908  FORMAT (IH  )  0066 

909  FORMAT(lHl)  0067 

FRELC 11 )s500.  0068 

FREL(12)«500*  0069 

DO  6461  IA*1»10  0Q70 

6461  I FR 1 1 A ) «FR€L ( I A )  007 1 

C  BEARINGS  OF  RMOMBICS 

IF(XETA)  1320*1320*1321  0072 

1321  XETAsABSFIXETA^XTR)  0073 

IFIXETA^IBO.)  1320.1320*211  0074 

211  XETAa360,*XETA  0075 

1320  IF(WiTA)  1323*1323*1324  0076 

1324  WETA«ABSF<WETA«BTRY)  0077 

IF(WETA»180*)  1323*1323*210  0078 
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2i0  WltA»36d,»WitA  6079 

1323  XQB.OOl  0080 

RQ>,00l  0081 

XP>4#  0082 

RP>4*  0083 

INH-XNH  0084 

INL«XNU  0083 

IND«XND  0086 

ONHaRNH  0087 

JNLsRNL  0088 

JNO-RNO  0089 

Ml 6x0  0090 

JiG-JlG^^l  0091 

HA>XA8SF(MAN)  0092 

C  CONSTANTS  TO  Bt  USED  FREQUENTUY 

TOP(l)>20i  0093 

TOP(  2)«20.  0o94 

TOP(3l*23i  0093 

T0P(4)«23.3  0096 

R0*6376«6  0097 

KIPal  0098 

IF(MaP)  36l2t36l2*56l3  0099 

3612  WRITE  OUTPUT  TAPE  6*3*  OlOO 

1  NOC I R  *  BMONS  *SSN  *  HA  *M AR  *  XL  A  T  *  AX  *  X  LONG  *  AL  »  YL  A  T  *  A Y  *  Y  LONG  *  YL  *  0 1 0 1 

2XTR*BTRY*GC0NM  0102 

3613  FX-69.037  0103 

tZP*!  0104 

LIP»1  0103 

MIP«2  0106 

IFIIHRI  8400*8400*3219  0107 

8400  GO  TO  (3219*3219*3220*3221«3219*3219) *NETH0P  0108 

3220  MIP>1  0109 

GO  TO  3219  Olio 

3221  MA>0  0111 

3219  DO  3214  |T«M!P*24*M1P  0112 

FREL(ll»*FOTnT)  0113 

JUGpO  0114 

KAT-l  0115 

DO  231  IJ«1*11  0116 

NANGLE(IU)»0  0117 


! 


4 
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MODI(IJ)>0  0118 

251  ABC(IJ)«AZZ  0119 

FRtQ>FREL( 1)  0120 

OFRfOi>3«  0121 

NF>1  0122 

C  F2  RlGiON  LAYER  HEIGHTS 

FL*0.  0123 

DO  926  II«XltK3  0124 

GALL  POLY  1 2 1 7  »  7  tGLCK ( I T )7 1 0. -1 *  2  »  GLAT  < 1 1 ) /90  / »  Y )  0125 

926  FL-Fl^Y  0126 

RX-KS  0127 

FLaFL/RX«100«  0128 

GO  TO  (927t928)*K2  0l29 

G  MODE  GhOOSING  F2  AND  E^LAYER 

928  ME(l)«GGDKM/2OO0.-f.9  0130 

ME(2)aME(li4>i  0131 

MF(l)«GGDKH/4000.-»i8  0132 

MF(2)*MF(U  +  1  0133 

MF(3)»MF(2)+1  0134 

GO  TO  929  0135 

927  ME<1»»1  0136 

ME(2)>2  0137 

MF(l»»l  0138 

MF(2)a2  0139 

MF(3)»3  Ol40 

G  Radiation  ancle  e-layer  modes 

929  DO  930  IK«lt2  0141 

1M>ME(IK)  0142 

P1>ME(|KI  0143 

ELDnM)«GGOXM/PI  0144 

EOELI  IM)iiATANF<  (GQSF  (GG0/(  2»4PI  )  I -RO/{RO-«-110»  ) ) /SINFI  GCD/ (  2«*PI  I )  )  0145 

1*BK  0146 

IF(EPiL(IM))  1032t93gt930  0147 

1032  DO  1033  II>lt2  0148 

1033  ME(|I )>ME( II )«1  0149 

GO  TO  929  0150 

930  continue  0151 

G  RADIATION  angle—  FZ-LAYER  MODES 

1934  DO  931  IK*1«3  0152 

1M>MF(|X)  P1S3 


-  63  - 


PI*iM  6iS4 

F  LD ( I M ) =  GCDKM/ P 1  0 1 55 

FOIL  < I M ) «ATANF ( <  COSF i 6CD/ ( 2 * *P1 ) ) -RO/ <  RO+FL  ) ) /S I NF ( GGD/ ( 2 . *P I » ) )  0156 

1*BK  0157 

IF(F0IU(IM))  932«931i931  0158 

932  DO  1933  11*1*3  0159 

1933  MF(  1 1  )=MFn  1  1  +  1  0160 

GO  TO  1934  0161 

931  GONTlNWl  0162 

GO  70(933*934) .K2  0163 

G  penetration  FREOUENGIES  -‘»E“REG10N 

934  0Z*EMF(4*1T)  0164 

0v*EMF(5tlT)  0165 

IFIOZ^OY)  940*940«941  0166 

940  ZRaOZ  0167 

12*2  0168 

2S*0y  0169 

GO  TO  939  0170 

941  ZRsOY  0171 

12*3  0172 

ZS*PZ  0171 

GO  TO  939  0174 

933  ZRsEmF(1*1T)  pits 

ZSsZR  0176 

939  DO  1663  1I=1»10  0177 

ELL(I1)>0.  0178 

1663  FLL( 1 1 )»0,  0179 

DO  827  11x1*2  0180 

IMxMEdl)  0181 

ERCR(IM)<Q,  0182 

EXMT<IM)s0.  0183 

<:  MODE  ELIMINATION  E  ANO  F2  REOIQNS 

ARZxELD( IM)*PK  0184 

1F(16*=ARZ)  2900*2900,2901  0185 

2900  ELF|i;<:l,02  0186 

GO  TO  2902  0187 

2901  ELFK=( ( « I ( (-4.368460907E-9*ARZ+1.33449426li”7)*ARZ-5,9T66l8436E"6)  0188 

l»ARZ+2,624e083i5E’'4)»ARZ-5, 0384762661-3 )*ARZ*3,761385Q53E*2>*ARZ=1  0189 

2, 1332007561-2 >»ARZ+, 2085  0190 

2902  Y=ZR*ELFIC+.05  0191 
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IFCY^^FREQ)  328«326«827 
328  ELK IM)- 1000^ 

827  CONTINUE 

DO  227  Ii>lt3 
IMaMF (11) 

FRCR ( IMI-O. 

FXMT( IH)>0. 

ARZ>RO*2.*(Pl2‘«FDEL(  lM)»AK'‘ASiNF(RO*SINF(Pl2-^FOEL(  lM)*AK)/(RO‘^110i 
1)1  )*OIC 

iFde.^^ARZ)  2903*2903*2904 

2903  ELFK«1*02 
GO  TO  2905 

2904  ELFKs( ( ( ( ( (-4*368460907E-9*ARZ-*i*33449426lE“7)*ARZ“5*9786i8436E*6) 
i*ARZ+2  *624808315E“4 ) *aRZ«5. 0384762662^3 ) •aRZ+S  *761 3850531*2 ) *ARZ*1 
2  *  133200756E*2  )  •aRZ-**  2085 

2905  y>ELFX*ZS4«05 
1F(Y»FREQ)  227*227.228 

228  FULl IMlAlOOO* 

227  CONTINUE 

00  946  Il<iK2*X3 

F24(  in»<F2S(lI.lT)»(l80*-S$N)+F2MUi  ♦iT)»(SSN*lO*)  )/l70. 

946  GMA( li )s(QHL( ll *iT)*( 180**SSN)-*OHH( 1 1 » I T ) *( $SN*10. ) )/l70* 
S1>(FRE0*GMA(X2) )/(F24(X2)*GmA(R2)) 

S2* ( FRE0*GMA ( K3 ) ) / ( F24 ( X3 ) *GMA ( X3 ) ) 

IF(S1-S2)  947*947.948 

947  R«$2 

GO  TO  949 

948  R*Sl 

949  DO  945  IK>1*3 
IM>MF( IX) 

TF(FLL(IM))  508*508*945 
508  ARX-FLD(IM)*DX 

FLF( IMI>(( ( ( ( (-6.712654756f*9*ARX+4,49l5144lE»7)*ARX*9,98583ll04E* 
16 )tARK*6. 8652598 17E-5)*ARX49. 2024373328*5 )<^ARX+2*264634341l«3)tARX 
2+4.6992431011=3 ) *ARX 
IF(FLF( )M)*R)  950*945*945 

950  FLL( IM)»1900* 

945  CONTINUE 

GYRri(GY(X2)+GV(X3)  )/2. 

GO  TO  (1740*9119) *X!P 


0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 
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G  MAJOR  LAND  BODIES 

1740  W'LD^O*  0232 

DO  826  n»<iiK3  0233 

IF  (GLONGdlll  1702#1702»1701  0234 

1701  CZ6s3604-CLONG( M >  0235 

GO  TO  1703  0236 

1702  GZG*A8SF(CL0NG(  ID)  0237 

1703  CAUL  N0lsy<6iX2#G2G*Y)  0238 

826  WLD»WUD+y  0239 

WDswuD/RK  0240 

9119  IF<wD)  953.954*954  0241 

953  ERsSO.  0242 

S1Gma*5.  0243 

GO  to  955  0244 

954  ER>4«  0245 

SiGMA*.001  0246 

G  E'^LAYfR  MODES 

G  GROUND  LOSS  REFLEGTION  FACtORS  AND  LOSS 

955  DO  726  1I»1.2  0247 

I M»ME ( I i )  0248 

iF(ELL<  IM) )600t600.726  0249 

600  RS*IM*1  0250 

CALL  LOSStEOELdMI.ER.FREO.SlGMA.Y)  0251 

6E<1M)»AB5F(Y»RSI  0252 

RI*IM  Q253 

G  sky WAVE  ABSORPTION 

ESKY(IM)«615,5»RI4ABI ( IT»/IC0SF<ASINF(R0*SINF(PI2+EDEL«'1M)*AK)/{R0  0254 
l-dOO.) I )*(FREQ+1.12»GYR)**1,98)  O255 

RS*IM  023^ 

C  ANGLE  AT  IONOSPHERE 

PHEpRO*SINF(PI2+eOEL« IM)*AK)/IR0+110. )  0257 

C  FREE  SPACE  LOSS 

ZOR-  GCP/(2.*RS>  0258 

FRs ( 2 . *RS*S I NF ( ZOR ) ♦RO/PHf ) •EEK  Q259 

FSE(  IM)B36.58-f20.*LOG10F(FR)dQ**L0G10F(FREO)  0260 

726  eONTINUE  0261 

C  F2-LAYER  MODES 

C  GROUND  LOSS  REFLECTION  FACTORS  AND  LOSS 

DO  626  II«1.3 


=  66  - 


1F«FUL(IM))  60 1*66 it 626  6264 

601  6265 

CALL  LOSS(FOtLMH)tER.FRiQtSIQMA»Y)  0266 

GF ( I M ) *ABSF  t  Y*RS )  6267 

RI-IM  6268 

C  SKY  MAVE  ABSORPtlON 

;FSKY< 1M)«615.5*R1*AB1( I f ) / (COSF < ASINF «RO*SI NF ( P I2»FDEL( 1M»*AK)/<R6  0269 

U160*) )  »*(FREQ+lil2*6YR)**U9Bl  6270 

RSBIM  6271 

C  ANGLE  At  IONOSPHERE 

PME*R0«SlNiF(PI2+F0EL(  IMJ*AK>/(R0»FL»  6272 

C  FREE  SPACE  LOSS 

ZOR«  GCD/(24*RS)  6273 

FR* ( 2 i*RS*SI NF ( ZOR ) *RO/PhE ) •EEK  0274 

FSFt lM)»36.58»2Oi*LOGl0F(FR»+20»*LO6lOF(FREO)  0275 

626  continue  6276 

C  ANtENNA  determination 

KAP«XABSF( IANT)  0277 

IF(IANT)  18  #181*181  0278 

181  KAP*4  6279 

DO  82  11*1 *2  6280 

IM«ME(1I)  6281 

82  iXMTl IH)>IANT  0282 

0062  il«l#3  0283 

IM«MF(M»  0284 

62  FXMt(lM)>IANT  0285 

GO  TO  1185  6286 

18  GO  TO  (31t32*33) tKAP  0287 

31  TANT(1)-ARA( 1)  0288 

TANT(2IsARA(21  0289 

60  TO  180  0290 

32  TANt(ll«AVA(l>  0291 

TANT(2I>AVA(2>  0292 

GO  TO  180  0293 

33  TANT(1)>AHA(1)  0294 

TANT(2)-AHA(2)  0295 

180  DO  182  ll8l*2  0296 

IM>ME(in  0297 

IFIELWIIM) )663«6Q3tl82  0298 

603  CALL  GAIN(KAP*EDEL(IM) .XETA*FREQ*XO»XP»XNDtXNL»XNHtEXNT(IMn  0299 
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182 

CONTINUE 

0300 

t 

POWER  GAIN  Of  TRANSMltTiNG  ANTENNA 

DO  382  II>lt3 

0301 

1M>MF( 1 1 1 

0302 

1 F 1  Flu ( i M ) )  605  *  605 • 382 

0303 

60S 

CALL  GAiN(KAP«FDEL( IM) tXETA*FREO*XQtXP*XND*XNL •XNH*FXMT< IM) I 

0304 

382 

CONTI NUE 

0305 

1185 

IF(METM0d«4)  185 t 184*185 

0306 

184 

IF(IHR)  185*185*183 

0307 

183 

lANRaO 

0308 

MA«0 

0309 

185 

KUP-XABSF( IANR) 

0310 

IF  (IANR)  187*281*281 

0311 

281 

KUP>4 

0312 

DO  482  It«l*3 

0313 

IM«MF(iI) 

0314 

482 

FRCR( iH)>!ANR 

0315 

DO  582  II>1*2 

0316 

IH>ME ( i I ) 

0317 

582 

|RCR( tM)«iANR 

0318 

GO  TO  1285 

0319 

187 

GO  TO  (741*42*43)  *KUP 

0320 

741 

RANT(1)>iARA(1) 

0321 

RANT(2)>>ARA(2) 

0322 

GO  TO  280 

0323 

42 

RANTd  )«AVA(  1) 

0324 

RANT(2)«AVA(2) 

0325 

GO  TO  280 

0326 

43 

RANT(1)«AHA(1) 

0327 

RANT ( 2 ) SANA ( 2 ) 

0328 

C 

RESPONSE  OF  RECEIVING  ANTENNA 

280 

DO  282  ni>l*3 

0329 

IMsMFdl) 

0330 

1F(FLL(|M)I  607*607*282 

0331 

607 

CALL  GAIN(KUP*FDEL( IN) *WETA*FREQ*RQ*RP*RND*RNL»RNH*FRCR( IK) ) 

0332 

282 

CONTINUE 

0333 

DO  782  II>1*2 

0334 

IM-MEdI) 

0335 

IF(ELL( IK) )6Q9*609»782 

0336 

609 

CALL  GAIN(KUP*EPEL( IM) *WETA*FREQ*RQ*RP*RND*RNL *RNH*ERCR ( IM) ) 

0337 

4 


782  CONTINUE  0338 

1285  IFIMAP)  285 *285*388  0339 

285  GO  TO  (286*388) *uiP  0340 

286  GO  TO  (330*330*330*1331) *KAP  0341 

1331  GO  TO  (232 *232*232*233 )*KUP  03A2 

233  wPiTE  Output  tape  6*i3*iant*ma*ianr  0343 

GO  TO  388  0344 

330  GO  TO  (335*335«335*336)*KyP  0345 

335  WRITE  OUTPUT  TAPE  6*1*  0346 

1  TANT«INH*INL.IND«HA*RANT»JNH*JNL*JN0  0347 

GO  TO  388  0348 

336  WRITE  OUTPUT  Tape  6*2*  0349 

1  TANT*lNM*lNL*jNO*MA*lANR  0350 

GO  TO  388  0351 

232  WRITE  Output  Tape  6*4*iant*ma*Rant*jnh*jnl*jnd  0352 

388  LlP>2  0353 

DO  333  II>iK2*<3  0354 

TIP«0.  0355 

C  OfTERMINATION  OF  TYPE  OF  ClRCUlT(POtAR*TEMPERAT|*fTC» 

RM«AeSF(GLAT(ItM  0356 

iF(RM’-60.*AX)  331*331*332  0357 

331  ID>1  0358 

IF(A6IY(i*lT))  9732*9732*333  0359 

9732  TIP*5.  6360 

GO  TO  333  0361 

332  IF(RM-70.*AK)  334*334.  35  0362 

334  1P>3  0363 

IF(K2-'l)  333*333*  36  0364 

36  10*4  0365 

GO  TO  1725  0366 

35  ID*2  0367 

333  CONTINUE  0368 

C  TRANSMISSION  LOSS  E^LAYER  MODES 

1725  DO  960  11-1*2  0369 

IM-ME(I!)  0370 

960  ELL(IM)«ESKY(  IM)<fGE(  IM)<fFSE(  IM|-<>EXMT(  IM)^ERCR(  IM)4APJ(  1P)4.|LL(  IMI  0371 

LD-ME(l)  0372 

LE-ME(2)  0373 

IF(ELL(LP)«ELL(LE) )  961*961*962  0374 

961  ELOS-ELL(LP)  0375 
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MDE*Ub  0376 

60  TO  963  0377 

962  ELOSsEUU(LE)  0378 

:M0E*WE  0379 

C  TRANSMISSION  LOSS  F2-LAyER  MODES 

963  DO  964  IIai»3  0380 

IM»MF(M)  0381 

964  FLL( IM)=FSKY( IM)+GF( IMj+FSFC IM) ^FXMf ( IM>*FRCR« lM)+ADJ(iD)+FLL( IM)  0382 

LG*MF(1)  0383 

LH*mF ( 2 )  0  384 

LJsMF(3)  0385 

iF<FUL<L6)=FUL(LH) )  96 St 965 >966  0386 

965  FLlaFLL(LG»  0387 

M0F=lG  0388 

GO  TO  967  0389 

966  FUiaFLL(LH)  039o 

MdF»LH  0391 

967  IF(FLi«FLL(LJn  968t968t969  0392 

968  FlOS»Fl1  0393 

GO  TO  970  0394 

969  FLOSaFLL(LJ)  0395 

M0F«LJ  0396 

970  IF(ElOS-FLOS)  97lt97lt972  0397 

971  NXL0$(NFI»EL0S  0398 

N ANGLE  (  NF )  •EpEL  ( MOE )  5  0399 

MODE ( NF ) «MPE  04OO 

ABC(NF)«AEC  0401 

GO  TO  911  0402 

972  NXLOS(NF»«FLOS  0403 

mangle (NF)sFDEL(MPF)+. 5  0404 

MOPE ( NF ) >MPF  0A05 

ABC(NF)*AFC  0406 

C  MIXED  MOPES  OF  PROFAGATION 

911  FHOP»MF(l)=l  0407 

IF(FHOP)  3903t3903t3e98  0408 

3898  XHOPaFHOP'flf  0409 

C  AVERAGE  height  OF  REFLECTING  LAYERS 

EFHT ■ ( 1 ig,*FM0P*FL ) / ( FHOP+l • )  OAlO 

C  RADIATION  ANGLE 


EFPEL»ATaNF( (COSF(GCP/(2.*XH0Pn^R0/tR0+EFHT) )/SINF(Gep/(2t4XHOP» » 


0911 


i)«BK  0412 

C  ANGLE  At  IONOSPHERE 

PHi»R0*SiNiF(Pl2+tF0lL*AK)  /(RO+ilO.  »  0413 

ARZs79ii38*(Pi2“EFDEL*A<“ASiNF<PHEn  04l4 

iFIEFDfL-l*)  3899 •3900. 3 901  0415 

3899  FHOP»iFHOP+l.  0416 

GO  TO  3898  0417 

C  MODE  EL I  M I NAT  I  ON 

3900  ELFK=i,02  0418 

Go  TO  3902  0419 

3901  ELFK»( I ( ( ( (^4.368460907E-9*ARZ+l,33449426iE"7)»AR2-‘5*9766l8436E“6)  0420 

i*AR2+2 *62480831 5E“4)*aRZ-5.G38476266E“3)*AR2i-3. 761 38505 3E^2)*ARZ“1  0421 

2*i3320O756i^2)»ARZ*-*20d5  0422 

3902  Y=ZS*ELFIC  0423 

4332  IF(V-FREO)  3903 .3905 *3905  0424 

3903  NOS(NF)«1000.  0425 

GO  TO  612  0426 

3905  YsZR*ElFk  0427 

IFIY-FREQ)  650*650*651  0428 

651  NOS(NF)i>1000*  0429 

60  TO  612  0430 

650  ARX«(6CDKM*0K-ARZ)/Fh0P  0431 

FLF<  #( ( < ( ( (^6.712654756E-9*ARlC+4*4915l44lE*7)*ARX-9*9e5831l04E»  0432 

l6)»ARK+6»8652598l7E-5) *49^+9. 202437332E-§)»ARh+2.264634341E-3)*ARK  0433 

2+4.6992431Q1E»3»*ARk  0434 

REF«(FRfO^QMA( 12  » l/( F2A( I2)-GMA( J2) )  Q435 

IF(FLFK,^REF )  652*653.653  0436 

652  NOS(NF >>1000.  0437 

GO  TO  612  0436 

C  SKY  WAVE  ABSORPTION 

653  EFSK Y =6 1 5 • 5*XHOP* AB I Y ( I . I T | / ( ePSF « AS I  NF  I RO*Sl NF I P 1 2 *EFDE L» AK ) / ( RO  0439 

1*100.) ) )*<FREQ+1,12*GYR)**1. 98)  0440 

C  GROUND  REFLECTION  LOSS 

3908  CALL  L0SS(EF0EL*ER. FREQ. SIGMA. Y)  0441 

3909  GEFsABSF(Y*FHOP)  0442 

C  ANGLE  AT  IONOSPHERE 

SPHE»RO»S! NF ( P I 2*EFPEL»AK ) / ( RQ+EFHT )  0443 

C  FREE  SPACE  LOSS 

Z0Rs6CD/(2.#(FH0P+1, ) )  0444 

FR« ( 2.« ( FH0P+1 , ) tSlNF ( ZOR )*R0/SPHE ) •££<  0445 
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E  F S F  s  5  6  «  5  8 + 2  0  *  * U d G 1  OF ( F  R ) +20 • *U0G 1 Of ( P  R  6  Q )  0446 

GO  TO  (3l80t3i80*3l80,3iei)»KAP  0447 

3181  EFXMf=lANT  0448 

GO  TO  3182  0449 

G  POWER  GAIN  Of  TRANSMITTING  ANTENNA. 

3180  GAUL  GAlN(  KAP«EFDEL»XiETA»F«EO»XQ»XP*XNDtXNiLiXNH»EFxMT  )  0450 

3182  GO  TO  (3183f3183«3183*3184)*XUP  0451 

3184  EFRGR=IANR  0452 

GO  TO  3190  0453 

G  RESPONSE  OF  REGElVlNG  ANTENNA 

3 183  Gall  ga i n  <  kup  *  eFoeu  *  wet  a  »  freo  tRo »  rp ♦ R  no  + rnu  t  rnh  *  eF  rgr )  0454 

G  Transmission  loss 

3190  N0S(Nf jsEFSKY+GEF+IFSF^EFXMT-EFKGR+AOJI 10)  0455 

612  IF(N0S(NF)=NXLOS(NiF)  )  623.622*622  0456 

623  NANGLE(NF)*EF0EL+.5  0457 

MOOE(NF)=FtNOP+l,  0458 

B  ABG(NF) 4676060606060  0459 

NXLOS ( NF ) =NOS ( NF )  0460 

622  1F(NXLOS(NF)-990  )  6911.624*624  0461 

624  ABG(NF)3A22  0462 

NANGLE«NF)=0  0463 

MOOE(NF)sU  0464 

GO  TO  8740  0465 

6911  KAJ4l  0466 

GOTO  (6912*6912*6912*9779*6912.6912) .method  0467 

6912  POR410.*LOG10F(PwR*1000.)  0468 

IFliuG) 3869*3869.8712  0469 

C  IMC/S  ATMOSPHERIG  NQISE  DETERMINATION 
3869  Ce=9KC(IT)  0470 

e  SEASON  AND  HOUR  BLOCK  FOR  imtOSKEERie  NOISE  DBrEIGgNATION  0471 

L!B=4  0472 

ICCaCC/2.  0473 

IF(ICC)  8889.8889*8888  0474 

8888  GO  TO  (8701*8701. 8702*8702. 8705*8705. 8706. 8706*8708*8708*8887.  0475 

18887). ICC  0476 

8887  CC»CC’-24.  0477 

8889  XJ4l  0478 

TM=-2.  0479 

JK^l  0480 

GO  TO  8703  0481 
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8701  KJ»i 
TMa2* 

JKif 

GO  to  8703 

8702  Kd»5 
tMs6. 

JK> 

KAJ 
60  TO  8703 

8705  KJ«4 
tiM>10« 

JKs3 

LIB«5 

60  TO  8703 

8706  XJ>3 

TMaiA. 

JK>2 

L1B*5 

60  TO  8703 
8708  ICJ>2 
Tm«18« 

JK>1 

KAJ*2 

8703  IF(Y2»  9702»9702»9701 

9701  GEG«360.-y2 
GO  TO  9703 

9702  CEG»ABSF<Y2) 

9703  CALL  N0ISY<KJtX2  .CEGtATNOJ 
CALL  N0iSY( JKtX2»GEG»ATNY) 

ATNQs  ( ATNO-*-(  ATNY^ATNO )  *  <  CC*TM )  /4,  » 

L0B«0 

JUG«JUG4-1 

GO  T0( 8712 *8713) tXAJ 

c  frequency  dependence  of  atmospheric  nqisi 

8713  M0T»M0USE/3 

IF(MOt)  8714*8715*8714 

8715  M0T»4 

8714  1F«X2)  8716*8717*8717 

8716  L0B>1 

8717  IF(MOT’-3)  8718*8719*8719 


0482 

0483 

0484 

0485 

0486 

6487 

0488 

0489 

0490 

0491 

0492 

0493 

0494 

0495 

0496 

0497 

0498 

0499 

0500 

0501 

0502 

0000 

0503 

0504 

0505 

0506 

0507 

0508 

0509 

0510 

0511 

0512 

0513 


0514 

0515 

0516 

0517 

0518 

0519 
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8718  UIB=5“L0B  0520 

SO  TO  8712  0521 

8719  LlBa4+L0B  0522 

8712  IF(UIB^4>  2179*2179,2178  0923 

2179  IF(ATNQ“20, )  2130*2178*2178  0924 

2130  ATNQ*20*  0929 

2178  GA4U  POLY(4iB*7*7*ATNQyiO*-9. ♦FREQ/10. ‘‘1*39*7)  0926 

ATNOsYAlO.+lSO.  0927 

G  galagtig  noise 

IF(FRB0“6MA( 2) )  8720*8720*872i  0928 

8721  6N0S=  K <(< l (( ,iO99O32130E*l0*FR£Q»*2442892795E-8)*FREQ-. 19854096  0929 

1081^^5 )  *FRE04‘,  1  513740943E“3  )*FRi0‘-.6306642l89E-2  )*FREQ+.i390178395E  0930 

2+0 )  *FREQ-.  1 70l088799i+l  )»FRE0‘*‘*  1121  i30396E»2)*FREO**3426399698i+2)  0531 

3*FRtQ+.2O20289494E+3  0532 

GO  TO  9797  0933 

8720  GN0S»lO0O.  0534 

G  MAN  MAOE  NOISE 

9797  IF(MAN)  8722*8723*8723  0939 

8722  KJ»XABSF(MAN)+9  0536 

S0B=O.  0537 

GALL  POLY  (XJ*i»10»$OB*FREQ*Y>  0538 

XNOIS»Y  0539 

GO  TO  8724  0540 

8723  XNOiS»MAN  0541 

C  OETERMINATION  OF  CONTROLLING  NOISE 

8724  RGNSE*ATN0  0542 

IF<ATN0*GN0SI  8729*8725*8726  0943 

8726  RGNSE=GNOS  0944 

8729  IF(RCNSE=XN0ISI  8729*8729*8728  0545 

8728  RGNSE^XNOIS  0546 

G  AVAILABLE  SIGNAL-TO-^NOISE  RATIO 

8729  XLOS*NXLOSINF)  0547 

ROT  ■  XLOS-POR  0948 

GOTsRCNSE^ROT  0549 

GO  TO(509*905*509*505»505*906)»METHOD  0990 

506  NXL05(NF)-GOT  0591 

GO  TO  9777  0552 

505  WANT»GOT*RSN  0953 

C  CIRCUIT  RELIABILITY 

IF(ABSF<WANT)«(T0P( I0)+TIP))867*867*869  0554 
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86#  Ri L  «  S I CNF (100.*  WAN  T )  0##  8 

GO  to  869  0556 

867  KJ^iO-t^iO  0557 

GO  to  (888t8733*8733*8733)*iO  0558 

88#  lF(ABIY(i.lt))  8733*8733.8732  0559 

#732  0560 

8733  GALL  ROLY(KJ*5*10*fREO/iO*‘‘1.35*WANt/iO**Y)  0561 

RIL>Y*10**-50*  0562 

869  lF<i6-i)  8734*8734.#735  0563 

8735  XRL390*  0564 

GO  to  8736  0565 

8734  XRls99.  0566 

8736  IF(REL-xRL)  6739*8739*873#  0567 

8738  REL>xRL  056# 

8739  IF(REL‘1*  )  8740*8740*8741  0569 

8740  RELaO.  0570 

8741  NXL0S(NF)«REL  0571 

KIP>2  0572 

IF(MEtH00-^3)  9777*3207*9777  0573 

G  ItERAtlON  FOR  LUF 

3207  IF(REL^90.)  3201 *3202 *3203  0574 

3101  GO  to  (3204* 3202 )*KAt  0575 

3202  FF(1T)«FRE0  0576 

GO  to  3214  0577 

3203  lF(FREQ«3f)  3208*3208*3209  057# 

3208  FF(it)?-3.  0579 

GO  to  3214  0580 

3209  FREQ«FRE0^.1*PFREQ  0581 

I F ( FREQ°3. ) 3208 •4850.4850  0582 

4850  KAt«2  0583 

GO  to  939  0584 

3204  ZJ«KJ  0585 

GO  TO  9777  0596 

9779  !F(IHR)  9777.9777*977#  0587 

C  FIELD  STRENGtM 

977#  ALPSSs«NXL05(NF)  0588 

NXLOS(NF)»ALOSS+107.2+8.6#594LO(i(?(FREQ)*.A,343*LOGF(lOOO.*PWR»  0589 

C  FURTHER  ITERATION  FOR  LUF 

9777  NFsNF+1  qS90 

IF(NF’<11)  914*914*260  Q591 


‘  is  ^ 

914  FRE6*fRfL(NF)  0392 

if (FREO-fRiL<ll > »  9l7,9i7»916  0393 

917  if (FREQ-30.)  939*939.916  0594 

916  NXLOS(NF)*0  0595 

FF(if)»0.  0596 

(30  to  9777  0597 

260  60  TO  I250t250*32l4*3212*7007t230> .METHOD  0596 

7007  GO  TO  ( 7009. 7008 ).tZ#  0599 

7009  WRITE  Output  tape  6.34.PwR.iRSN.IFR  0600 

7008  WRITE  OUTPUT  TAPE  6.94.  0601 

1  iG( If ).uF({T).F0T(iT).NXUOS  0602 

Nl6sNlG-.^l  0603 

GOTO  (3214. 3214. 21) *NlG  0604 

2i  WRITE  Output  Tape  6.908  o605 

NIG>0  0606 

GO  TO  3214  0607 

250  GO  TO  (2250*2666) .LZP  0608 

2250  WRITE  OUTPUT  TAPE  6.14.PWR*IRSN.IFR  0609 

2666  GO  TO  (2667.2667.2667*2667.2667*2668) .METHOD  0610 

2668  WRITE  OUTPUT  TAPE  6*29.  0611 

1  IG( lT-l)*UF(lT-l).FOf(IT«l)*(HODE( 1 ) *ABC( I ) .I«l *11 ) .MANGLE  0612 

2 . IG( I T ) *UF (IT) *FOT ( 1 T ) .NXLOS  0613 

GO  TO  3216  0614 

2667  WRITE  OUTPUT  TAPE  6.9.  0615 

1  IG( I T*l) .UF( IT-1) *FOT( IT-l) .(M0DE( I).ABC( I) *I«l.ll).NANGLE  0616 

2 . IG( I T ) .UF (IT) .FOT ( I T ) .NXLOS  0617 

GO  TO  3214  0618 

3212  GO  TO  ( 3213*3215) *LZP  0619 

3213  IF(IHR)  8002 *8002 '8003  0620 

8003  WRITE  OUTPUT  TAPE  6»8Q*PwR*IFR  0621 

GO  TO  77  0622 

8002  WRITE  OUTPUT  TAPE  6*8* |FR  0623 

3215  IF(IHR)  69*69*77  0624 

69  WRITE  OUTPUT  TAPE  6*19.  0625 

1  IG(IT»1)*UF( IT-1)*FOT(|T>1)*(MOOE( I). ABC(1).I«1*11). mangle  0626 

2.IG( IT)*UF( IT) *F0T(IT)*NXL0S  0627 

GO  TO  3214  0628 

77  WRITE  OUTPUT  TAPE  6.59.  0629 

1  IG  ( I  T*1 )  .UF  (IT-1)  *FOT  ( I T-1 )  » ( MODE  <  I )  .ABC  (D.Ul.U).  MANGLE  0630 

2 . IG ( I T ) «UF ( I T ) »FOT (IT) .NXLOS  0631 
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3214  LZP>2  0632 

GO  to  (2270«31*290t51i41t31)tMttHOD  0633 

41  GO  TO  (40*31)* JIG  0634 

40  WRITE  OUTPUT  TAPE  6*907  0635 

GO  to  270  0636 

31  WRITE  OUTPUT  TAPE  6*909  0637 

JIG-0  0638 

GO  TO  270  0639 

290  IP (MAP)  6228*6228*6229  0640 

6229  CALL  CURVY ( POT •PP)  0641 

GO  TO  270  0642 

6228  WRITE  OUTPUT  TaPE  6*24«PwR*IRSN  0643 

IP(fHR)  2291*2291*2290  0644 

2291  WRITE  OUTPUT  TAPE  6*7.  0643 

1  (lG(I)*PP(I)*F0T(n*lG(I-*6)*PP(I4^6)*P0T(I-t-6)*IG(l4l2)*  0646 

2  PP(  H12)*P0T(  I-.UI.IGI  i*-l8)*PP(t-*18)  *POT(  I*-18>*  1-1*6 )*AAA  0647 

GO  TO  2270  0648 

2290  WRITE  OUTPUT  TAPE  6.70*  0649 

1  ( IG(  I)  *P0T(  1  )*IG(I4^6)*P0T(  I-*6)*IG(  I-*12)*P0T(  I-*12)*lG(  I^-l  0650 

28)*FOT  (  1+18)  *I6(  1-i)  *PP(  !♦!)*  POT(  I  +  l )  *I6(  I<-7)  *PP(  1+7 )  *POT(  U7>  *IG  0651 

3(  I+13)*PP(i  +  13)  *POT(  I*-i3  )  *16(  I  +  IOlPF)  1  +  19)  *POT(  1  +  19)  *1-1  *6*2)  *AAA  0652 

2270  IP(3-MIT)  1969*1969*1970  0633 

1969  WRITE  OUTPUT  TAPE  6*909  0634 

MIThI  0633 

GO  TO  270  0636 

1970  MlT«HIT+i  0657 

270  RETURN  0638 

END  0639 
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SUBROUT  INI  CURVY ( COT • BOf )  000 1 

C  GENERATES  UINE  GRAPHS 

OIMENSION  C0T(24) tROT(73) tBOT(24) *WOt(73)  0002 

DIMENSION  AOT(73) rCOB ( SB ) *B08( SS ) *XOB ( 73 )  0003 

OIMENSION  FOEKdOliEDEKdO)  0004 

DIMENSION  AM0N(l2)t  GLAT ( S ) »AB1 ( 24 ) • ABl Y ( 3 >24 ) #CUCK ( 24 ) tGY ( 5 ) *  0003 

lf2S(3t24) if2H(3t24) tGMU(3>2A ) tGMH ( 3 t 24 ) tMF ( 3 ) >ME(2 ) >EUL( 10) > FLU (10  0006 

2  )  >  E  U  D  d  V  )  )  E  D  E  k  (  1 0  )  >  F  L  D  d  0  )  >  FOE  L  ( 1 U  )  >  F  2  4  (  3  )  .  G  M  A  (  3  )  t  F  L  F  (  1 0  )  >  F  SK.  Y  d  0  )  0 00 7 

3>£SICYdO)  tGEdO)  >GF  (  lU)  >FSEdO)  >Fi>F(  10  )  tNOSdl  )  >FXMt  (  10  )  >EXMT  (  10)  >  0008 

4FRCR ( 10 ) >ERCR ( 10 ) t ADJ ( 4 ) > ABC ( 12 ) >MOUE ( 12 )  0009 

DIMENSION  Q(20>60>4)  tiL<4)  >J'U(4)  >lkL(4  )  >UX(4)  >JAL(4)  •CXC(24)  0010 

ltG(60) iAB(60) tS(20>24) tC(20>24) >GAMMA(3>24) >BA(60)  0011 

DIMENSION  R0(5) tSUN( 12 ) >A ( 1C >7 > l4 ) >RASSN (12 ) > EmF ( 3 >24 )  0012 

dimension  CLAT(3) >CU0NG(3)>EMUFY(3>24) >FMUFY(5>24) >  UFY ( 5 ) >F0T Y ( 3 )  0013 

3>IG(24)  0014 

DIMENSION  UF(24)>F0T(24)>FF(24)>P(29>16>&)>ABP(2)6)  0013 

DIMENSION  AHA(2)>ARA(2)>AVA(2)>RANt(2)>TANT(2)  0016 

DIMENSI0NNXU0S(  11  )  >NANGLEdl )  0017 

DIMENSION  T0P(4) >FREU( 12)  0018 

COMMON  ABP  >P  00 1 9 

COMMON  A > U  > G  >  AB >  S  > C  > GAMMA  >  BA  >  EMUF Y  >  F MUF  Y  > I G • U F  t  Fu T  > SUN  >  A V  A  >  AR A  > GL A  0020 

lt>ABI >ABI Y>CLCX>GY>F2S>.F2H>GML>GMH>MF*ME>IL0>FlP>F24>QMA>ADj>ABC>  0021 

2EMF>CKC>CLAT>CLONG>AX>B<>CK>OK>PI2>EK>EEK.>GLt tGLG>AAA>ALA>ASA>  0022 

3AFC>ANC>AWC>AEC>ASC>A2Z>MAN>NSN>KM>PWK>IANT>IANK>y2>X2iM0USE>SSN>  0023 

4NGCrR>X4>  AX>XLONG>At>Yi^AT>XLAT>YLONG>YLtBTRY>GCUXM>GCPNM>M>X2>K  0C24 

5 3  > K 5  > GC 0 > I R SN > MON > I p  >MOUSE >AY>METMPP, AHA > BMONS  » MA P  0023 

COMMON  FREU  0026 

COMMON  XNH>XNI^>XNP>RNH>RNL>R:NO>HA>U!G>HAR>XTR>XETA>WETa>IHR>MIT  0027 

COMMON  XAP>KUP>MA>TANT>  INh>INL>  INUtRANT>  JNH>J.NL>JNP  0028 

73  format (lOX>4HANT=13)2HDe>2CX>6HNpiSE= 13 >20X>4HANT=I3>2HOB)  0029 

91  F0RMAT(6X>2A6>1X> I3>1HH>1X>I3>1MI.>1X>  t3>3HpEG>2Xi6HNClSE=13>2A6>lX  0030 

1 > ! 3  > 1MH> IX  > I  3  > IHL  > IX  > 13  >3HPEG)  0031 

72  F0RMAT(6X>2A6>1X> I3>1HH>IX> I3>lHL>lX>I3>3HPiG>12X>6HNOtSI’|3>12X>  0032 

14HANT»I3>2HDB)  0033 

74  FORMAT! 10X>4HANT=13>2HDB>13X>6HNOISE=I3>7X>2A6>lX>I3>lHH>lX>I3>iHL  0034 

1>1X>13>3HPEG)  0033 

24  FORMAT! igx>4HPWR=F6,2>2MXW>42X>8HRIU,S/N=I3>2HOB)  0036 

3  format (14X> I S>10X>A6>  10X>4HSSN=F5,Q>10X>A2 >F7.3/13X>  IIHTR  0037 
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2  ANSM I  T  t  er  •  1 3 X  t  BHR ICE  I V »  1 2 X  » SHBE AK I NG S  •  6 X  »  7 HN  . M j  L E S / 1 1 X  .  f  6  •  2  »  A 1:  >  00 36 

32M  -  tF7.2*Aii5XtF6»2iAl»2H  -  »F 7. 2 . Al f3X » 2F6. 1 t4X t Ftt *  1)  0039 

15  FORMAT (7X»73Al/6Xt74HQO  02  04  06  06  10  12  14  0040 

1  16  18  20  22  00/34X«19MGREENWiGH  MEAN  TIME)  0041 

2  FORMAT(2XtAl«lXiA2f A1 i73Al*Al»A2«lXtAl )  0042 

25  FORMAT (IHlV  0043 

i  FORMATUH  )  0044 

WRITE  OUfRUf  TAPE  6.3*  0045 

1  N0CIRiBM0NS*SSNrHA»«AR»X4AT»AX.X40N6»AU»Yl.AT*AY»YLGN6»yLi  0046 

2XTR*BTRy*GCDNM  0047 

IFIMETM0D~3)  4044*4046.4045  0048 

4045  GO  TO  (330*330*330.1331) .KAP  0049 

1331  GO  to  (232*232«232*233) *XUP  0050 

233  wRlTE  OUTPUT  TaPE  6*73* I ANT.Ma* IANR  0051 

GO  to  388  0052 

330  GO  to  (335  •335*335*336  )*ICUP  0053 

335  WRITE  OUTPUT  TaPE  6.91*  0054 

1  TANT*INH*iNL*tNO*MA*KANT»JNH*JNL*JND  0055 

GO  TO  388  0056 

336  WRITE  OUTPUT  TaPE  6*72*  0O57 

1  TANT*INH*INL*IN0*MA*IANR  0058 

GO  TO  388  0059 

232  WRITE  Output  Tape  6*74*  0060 

1  IANT*MA*RANT*JNM*JN(,*JN0  0061 

388  LlP»2  0062 

WRITE  OUTPUT  TAPE  6*24*PWR*IRSN  0063 

4044  J»0  QQ^4 

00  303  1-1*24  00^5 

X-eOT(I-l)  0064, 

IF(I-I)  202*203*202  0067 

203  X-eOT(24)  0068 

R0T(l)sx  0069 

202  XX-eOT(I)  0070 

DO  303  11-1.3  0071 

0072 

CIO- II  0073 

303  ROT( J+l)-X+(XX^X)»CEG/3,  0074 

IF(METMOP^3)  6402*6405*6402  0075 

6405  J-0  0076 

DO  606  l9]Lf24  0077 


4 


-  79  - 


XsXX 

0078 

IF (1*1)  8u2iB03»8G2 

0079 

803 

X=®0T(24) 

0080 

tF(X)  5lt53«54 

0081 

51 

X  =  U 

0082 

60  TO  54 

0083 

53 

X=COT(24) 

0084 

54 

WOT ( 1 ) =  X 

0085 

802 

XX*BOT(i) 

0086 

IF (XX)  55 f 5 7, 58 

0087 

55 

XX>1. 

0088 

GO  TO  58 

0089 

57 

XX=OOT< I ) 

0090 

5a 

DO  808  iMlf3 

0091 

JsJ-f  1 

0092 

ge6» i 1 

0093 

808 

W0T(  J+1  )=X-t-(XX»X)  *CE6/3i 

0094 

6402 

DO  608  U2.72.1 

0095 

fi  60  8 

XOB( 1 ) >406060606060 

0096 

DO  20  I=1.73#3 

0097 

B  20 

XOB( 1 ) >206060606060 

0098 

DO  30  1>1«55 

0099 

8 

BOB( I )«606060606060 

0100 

o 

OD 

COB( I ) >606060606060 

0101 

B 

C06(l)>  6U0 3 60 60 6060 

0102 

B 

COB ( 5 ) -  600460606060 

0103 

8 

C0B(9)>  600560606060 

0i04 

B 

COB ( 13 ) >600660606060 

0105 

B 

COB( 17 ) >600760606060 

0106 

B 

COB ( 2 1 ) >60 1 060606060 

0107 

B 

COB( 23)>601160606060 

0108 

B 

COB ( 25 ) >01 1260606060 

0109 

B 

COB( 27 ) >010160606060 

QUO 

B 

COB( 29 ) >010260606060 

0111 

B 

COB( 31 ) >010360606060 

0112 

B 

COB ( 33 ) >010460606060 

0113 

B 

COB( 35 ) >010560606060 

0114 

B 

COB( 37 ) >010660606060 

0115 

B 

COB ( 39) >010760606060 

0116 

B 

COB ( 4 1 ) >  0 1 1 060606060 

0117 

B 

C&B( 43 ) >01 1 160606060 

0118 

B 

COB ( 45 ) >02 1 260606060 

0119 

B 

COB ( 47 ) >020260606060 

0120 

B 

C0B( 49 ) >020460606060 

0121 

B 

C0B( 51 ) >020660606060 

0122 

B 

COB (53) >02 1 060606060 

0123 

B 

COB ( 55 ) >031260606060 

0124 

B 

BOB ( 29 ) >456060606060 

0125 

B 

BOB (21) >622060606060 

0126 

B 

BOB( 22 ) >616060606060 

0127 

B 

BOB( 23 1 >236060606060 

0128 

B 

B0B( 24 ) >446060606060 

0129 

B 

BOB ( 27 ) >706060606060 

0130 

B 

BOB( 28 ) >236060606060 

0131 

B 

BOB( 30 ) >256060606060 

0132 

B 

BOB ( 31 ) >646060606060 

0133 

B 

BOB ( 32 ) >506060606060 

0134 

B 

BOB( 33 ) *256060606060 

0135 

B 

BOB( 34) >516060606060 

0136 

B 

BOB ( 35 ) «266060606060 

0137 

B 

OOf >546060606060 

0138 

20B>31« 

0139 

B 

OOB>406060606060 

0140 

DO  200  IK* It 55 

0141 

iM>56*IX 

0142 

DO  100  |>lt73 

0143 

8 

100  AOT( I ) >606060606060 

0144 

lF(Z0B^20t)  4t4t63 

0145 

63  N0B>1« 

0146 

QQ  TO  5 

0147 

4  )F(ZOB’>8t)  4062t4062t4063 

0148 

4063  N0B>.5 

0149 

60  TO  5 

0150 

4062  WOBs.25 

0151 

5  ZQB>ZOB«WOB 

0152 

IT»1 

0153 

12  ISPOT«rT 

0154 

|F(ROT(  ITj^^ZOB)  6t8t8 

0155 

8  SOB>ZOB-fWOB 

0156 

!F(ROT( !T)«SOB)  10f6t6 

0157 

-  81  * 


10 

AOT 1 i SPOT  I eOOT 

0158 

lf«IT+2 

01S4 

GO  TO  16 

0160 

6 

jf»IT‘*.l 

0161 

16 

lF(iSPOT“73)  i2»l3*l3 

0162 

13 

If  <MeTH00*3 )  200r6565»2Gy 

0163 

6565 

IT»1 

0164 

33 

ispot=it 

0165 

If (WOT( It)-ZOB)  66*38*38 

0166 

38 

SOBaZ0B4-WOB 

0167 

If (WOt(lT)^SOBl  80*66*66 

0168 

80 

If  (AOTI  ISPOT)‘‘OOT  )  81*66*81 

0169 

81 

AOT ( I  SPOT ) >336060606060 

0170 

ITalTfZ 

0171 

GO  TO  67 

0172 

66 

IT>IT*1 

0173 

67 

If<lSPOT-73l  33*200,200 

0174 

200 

WRITE  OUTPUT  TAPE  6*2* 

0175 

1  B0B( IM) *C0B( tM>*D0B*AOf *D0B*C0B( IM) *B0B( IM) 

0176 

14 

WRITE  OUTPUT  Tape  6*15*X0S 

0177 

WRITE  OUTPUT  Tape  6*25 

0178 

RETURN 

0179 

INP 

0180 
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SUBROUt I Nf  G A I N ( KOP  t  DE  L  T  A ♦BE  t  A ♦ FMC ♦ S I GMA ♦ ER ♦ PH  I ♦ EL ♦ M ♦ R A IN)  0001 

G  POWER  GAIN  OF  ANTENNAS 
97A9  format  (iFi2*6) 

RM I « PH I  * .01 74532 9252 

RELT A-OfLTA* ♦O 1 745329252  ®®®^ 

GO  TO  )6f7.6)tK0P  ®®®* 

6  T8C0SF(RELTA)  ®®®* 

RETA«BETA*.0i745329252  ®®®^ 

0»SlNF(RiLTA)  ®®®® 

R»Q»Q  ®®®^ 

S»R«R  ®®^® 

X*18000,#SIGMA/FMC  ®®^^ 

RH0«S0RTFI (ER-T*T)*(ER*t*T)»X*X)  ®®i2 

RH0l2aSQRTF(RH0)  ®®^^ 

ALPHA8-ATANF  <X/ 1  ER-^t*T ) )  ®®^^ 

PS 1 H* A  T ANF  <  2 . •RhO 12*0*  S I NF (ALPHA* . 5 ) / ( RhO^R ) !  001 5 

CH8SQRTF(RH0*RH0-*'S«-2.*RM0«R*C0SF|  ALPHA)  )/IRH0+R+2»*Rh012*Q*GGSF(AL  0016 

1PHA8.5I)  ®®^^ 

WAVE8299.7925/FMC  ®®^® 

FAG8 (3#l4l592654*EL) /WAVE  ®®1’ 

Ul*l.»T*SiNF(RHI*RETA)  0020 

U28l.*T*SlNF(RH|^REtA)  0021 

6FF80.  ®®22 

X*l.  00^3 

GO  TO  (2t2.31fKOP  ®02R 

2  R AI N- (3 , 2*C0SF ( RH I ) *COSF ( RHI ) *S I NF C F Ae*U 1 ) *S I NF ( FAG*U 1 ) *SI NF ( FAC*U  0025 

12 ) *S INF ( FAC*U2 ) * ( CH*CH+l.*2«*CH*C0SF ( PSIH'^  ( 12*  56637062*H) /WAvE*9) I  0026 
2  )^*01*U1*U2*U2)*(C0SF(RET4)*SINF(RMI )*T>*(G0SF(RETA)"^SINF(RHI  )*T)  0027 

R0K«3.  ®®?® 

60  TO  4  ®®^’ 

3  RAIN«(  •5562474*eOSF  ( 1,570796*1  )*COSFI  1 .570796*7  )*(CH*CH'*-l.»2.*CH*C  0030 

10SF(PSIH»( 12.56637062*H)/WAVE*Q) I )/S  0051 

R0K8*4,77  0032 

GO  TO  4  0053 

7  E»ER  0054 

IF(H)  24. 24,20  0055 

20  WAyE8299,7925/FMG  0056 

A>(6.285185*H)/WAVE  0057 

X8H/WAVE  0058 

60  TO  25  0059 
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24  EN»INTF(ABSf 0040 

IF(EN)  50#51t5O  0041 

51  M1«ABSF(H)  0042 

60  TO  52  0043 

50  H1«EN/< ABSF(HJ  »10.*iN>  6044 

52  AaHl«6.283165  0645 

X*Hi  0046 

25  0*2i*A  6047 

ROK«+5»  0648 

2«2i*0  6049 

W«0i(ZtWl)  6656 

W2«C1(0*W3)  0051 

RA»30.*I=*5»COSF<D)*(.5772156649+L06FI2)*W  )♦( l.+COSF ( Dl ) *( *577  0052 

12156649i-l06F(0|‘‘W2  >4S1NF<0)*(  •5»W1«W3  > )  0053 

26  S>SiNF(REL7A)  0054 

C*G05F(RELTA)  0055 

HM«A/6«283185307  0056 

ANyM»C0SF(A*SJ^COSF(A)  0057 

28  call  vRiFC0(RELTA«E«FMC*St6MA*CV»PSlV)  0058 

FAC1*CV*CV+1.-!‘2«»GV*cOsF(PS1V-12,5653706*MM*S>  0659 

FAC2»2**C6$F(6*2d31853072*MM*S)  0066 

29  RAIN  »( 120«*ANUM»ANUM)/<RA*C*C>  0061 

RAIN>RAiN«FACl  0062 

4  IF(RAIN)  11*11*31  0663 

31  RaIN«4.342944819*LOGF(RAIN)  0064 

IF(X^*25)  69*70*70  0065 

69  EFF»( ( «64l6.702573*X-6091.333295)*X+2l79,e90548l»X^364;8173803)*X+  0066 

125*64620146  0067 

GO  TO  71  0068 

70  EFF-0.  0069 

71  RAIN»RAIN»EFF<-ROX  0070 

41  IF(RA!N-*10.  )11«11*10  0071 

11  RAIN»-1Q*  0072 

10  RETURN  0073 

END  0074 


ee9""</ 


SUBROUTINE  VRlFCO ( OtU T A t ER tFREO tS IGMA  tCV «PS 1 V )  0001 

G  VERTICAL  ground  REfLEGtlON  GOtFflCiENT 

Piaa* 141592634  0002 

X  > 1 8  000  *  4S I GMA /FREQ  0003 

U» ( ER»ER+X*X )  0004 

VsSQRfF<U)  0005 

Q«StNF(OElkTA)  0006 

R«0*0  0007 

S«R»R  0008 

T*C0SF( DELTA)  0009 

RH0«SQRTF(  (ER-T*T)*(ER‘»T*T)+X»X)  0010 

,  RH012aS0RTF(RH0>  OOll 

ALPhA*-ATANF(X/<ER^T*T))  0012 

A*2.«RHOl2*Q«V»S!NF(ALPHA**54'ASiNFIX/V)l  0013 

B«RH0-U»R  0014 

iF(B)  10*20*30  0015 

20  if (A)  21*22*23  0016 

21  PSIV«*.5»P1  0017 

GO  To  50  0018 

22  PStViO.  0019 

GO  TO  50  0020 

23  PSiV«*5«Pl  0021 

GO  TO  50  0022 

10  P$1V«ATANF(A/B)4PI  0023 

GO  TO  50  0024 

30  PSiV-ATANF(A/B)  0025 

50  CVsSQRTF ( RH04RHO+U*U»S»2  «  #RHO*U*R*COSf ( ALPHA+2 . *45 1 NF U/ V  )  ) ) / ( RHO+  0026 
IU«R+2»*RM012«V*Q»COSF(ALPHA*.5*-AS1NFix/V) ) )  0027 

RETURN  0028 

END  0029 
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FUNCTION  CKX.Sl) 

C  SINE  AND  COSINE  INTEGRAL 

4  IF (X-I0i)5. 60*60 

5  S0»X*X 

Ci*.5772i56649+LOGF(X)  ®®®‘^ 

TERM»-i**SQ/4*  ®®®^ 

6*4*  0®®6 

10  Cl-CI+TERM  ®®®^ 

TERMs*1**TERM*SQ»(6-2*  »/  l  <G‘‘l*  »*Gt6)  0008 

Gs6+2*  ®®®^ 

IF(ABSFtTERM)“.00005)  20.20*10  ®®i® 

20  SI-0.  ®®“ 

TERM-X  ®®^^ 

6-3*  ®®13 

25  SI-SI+TERM  ®®^^ 

TERM»‘‘i**TERM*SQ*(6^2*  l/<  <6“l*»*6*Gj  ®®15 

6«6-t.2*  ®®^® 

IF(ABSF(TERH>“*00005)  80*80*25  0017 

60  TERMl-l*  ®0^* 


C-1.  ®®2® 

61  FAC-(4**T»T-2**T)/<X-X)  P02i 

TgRH2-^l.*FAC»TERMl  0022 

IF(ABSF(TERM2)-ABSF<TeRM1) 162*64*64  0023 

62  C-C+TERM2  0024 

T.T*1*  0025 

IF ( ABSF { TERM2/G  >  *10.E-10 » 65*65*63  0026 

63  TERM1-TERM2  0027 

GO  TO  61  0028 

64  SC-2**T  0029 

THETA-X*SC  0030 

PI»*5-< l.*2.*THeTA|/U,*5e)*«l*'4.*THETA-2**THETA*THiTA)/<8.*SC*SC  0031 

1)+(3.+18.»THETA+8*»Theta*THETA»/( 16.*SeASC*Se)  0032 

C«e+TERM2*PI  0033 

65  TERMl-l./X  0034 

P-TERMl  0035 

T-i.  0036 

66  FACs(4.*T*T+2*»T)/«X*X)  0037 

TERM2--l,»FAC*TERMl  0038 

! F <  ABSF ( TERM2 ) ^ABSF <  TERMl » J67 *69*69  0039 


$- 
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67  0>b-fTERM2  0040 

Tpf+li  0041 

IF  ( ABSF  ( TiRH2/D )  ‘^iO.EaiO )  7O*70*68  0042 

68  T|RMlPfERH2  0043 

60  TO  66  0044 

69  S0a2**T+l.  004S 

tHetA»X*SO  0046 

PI*,5‘"(li^2p»TMETA)/(4.*S0»»«i^*4#*THETA^2i*fHETA*THETA)/<8**SD*S0  0047 

l>4-(3.4-i8.*rHETA4^8.«tMETA*  THETA) /(16**SO*S0«S0)  0048 

Dp04.TERH2*PI  0049 

70  CX»€0$F(X)  0050 

SX>S1NF(XI  0051 

G|«(6«SX*D*CX)/X  0052 

SI«i*570796327^(C*GX4-0*SX)/X  0053 

80  return  0054 

END  0055 


! 

1 
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SUBROUTINE  L0SS(RELtA*ER*RREQtSi6MAtY I  0001 

GROUND  reflect  ION  LOSS 

fELtA«RiLtA«. 01745329  0002 

PI >3 i 141592654  0003 

XalaOOOi«SlGMA/FREO  0004 

U«  ( ER*ER‘«‘X*X )  0005 

V>SQRTF(UI  0006 

Q«SInF( DELTA)  0007 

R>GN»Q  0008 

S*R*R  0009 

T>COSF (DELTA)  0010 

RHOaSORTFl (ER-T*T)*(ER*T*T)+x*X)  0011 

RH012aSQRTF(RH0)  0012 

ALPHA**ATANF ( xy ( ER-T*T » )  0ol3 

CHaSORTF ( RH0»RH0+S^2 .4RH0*R*C0SF ( ALPHA ) ) / 1 RHP+R+2  a*RHOi 2*0*C0SF ( AL  00l4 
1PHA*.5))  0015 

50  CVaSQRTF(RMO»RHO+U*U*S^2.*RHO*U4R*cOSFtALPHA*2»*ASINF(X/V)))/lRHO+  0016 
1 U*Ra  2 • *RHO 1 2  « VAQftCOS  F ( ALPHA* • 5 ♦ AS I NF I X / V ) ) )  OO 1 7 

Y«l0.*LO6i0F( (CH*GH+CV»CV)/2*)  0016 

REtURN  0019 

END  0020 


J 


4 
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SUBROUTINE  POLY( KJ>NNtHH*V*X»Yt 

POWER  SERTES  VARIABLES 

oool 

DIMENSION  A(L0»7*iA) 

0002 

OIMENSION  P( 29»16*6) *ABP( 2*6) 

0003 

COMMON  ABP*P 

0004 

COMMON  A 

0006 

Y*0 

0006 

MiiMM 

0007 

120  N-NN 

0008 

COEF-O 

0009 

127  COEFaCOiF*V*A(M*N«KJ) 

0010 

N-N-1 

0011 

tF(N)  122*122*127 

0012 

122  Y«Y*X*C0EF 

0013 

M>M^1 

0014 

tF(M)  126*126*120 

0016 

126  RETURN 

0016 

END 

0017 

-  89  - 


SUBROUt INE  NOISY (KJ*XP*CiGtATNO) 

COURIER  VARIABLES  AND  AtMOSPHERlC  RADIO  NOISE 

0001 

OIMENSI ON  PI  29 « 16 •& ) • ABP ( 2  »6 ) *2Zl 29 ) 

0002 

COMMON  ABPtP 

0003 

AUF»ABP(itKJ) 

0004 

BEtaABP(2iKJ) 

0005 

0>0«008t2664664626*CE(> 

0006 

CiscOSFiO) 

0007 

S1>SINF(0) 

0008 

DO  56  J>1»29 

0009 

R»0« 

0010 

Sx>si 

0011 

CX«C1 

0012 

DO  55  K«l«i5 

0013 

R>R4-SX*P|J*K*XJ) 

0014 

SS«SX*C1-»<X«S1 

0015 

CX*GX*Ci-Sx*Sl 

0016 

sx>ss 

0017 

Z2<Jt4R4‘PlJ«i6»KJ) 

0018 

Q*«0l745329252»f  XP4-90. ) 

0019 

Si>SlNF(0) 

0020 

SX>S1 

0021 

Cl«C0SF(0) 

0022 

CX«C1 

0023 

R»0* 

0024 

DO  57  K>lt29 

0025 

R»R*SX»2Z<X) 

0026 

SS«SX*C1+CX»SI 

0027 

cx«cx*ci»sx*si 

0028 

sx«ss 

0029 

ATNO-R-»'Al.F-i«Er*0 

0030 

RETURN 

0031 

END 

0032 

IX.  MATHEMATICAL  EXPRESSIONS 

L  Oreat  circle  distance  in  kilometers,  degrees  and 
nautical  miles 

cos  (gcd)  -  sin  (xl)  •  sm  (x2)  cos  (xl) 
cos  (x2)  •  cos  (y2  «  yl) 
gcd  -  great  circle  distance  -  degrees 
xl  =  transmitter  latitude  -  degrees 
yl  =  transmitter  longitude  -  degrees 
x2  =  receiver  latitude  -  degrees 
y2  =  receiver  longitude  -  degrees 

2,  Bearing  from  receiver  to  transmitter,  degrees 
east  of  north.  (0  ^  btry  —  360*). 

btry  =  114.  5816  •  tan"^  •si .  ;  ,  _  ,  .  ..i  /  i  j 

sin  (u*90  4x2)  ‘  sm  (u«ged)  /cind 

u  s  (180.  *  xl  •  x2  4  gcd)/  2 
cind  -  sin  (u)  •  sin  (u  «  90  4  xl) 

3,  Geographic  latitude  of  control  points  along  great 
circle  route. 

x'  1  =  cos  (pp)  •  cos  (90  -  x2)  +  sin  (pp)  •  sin  (90  xZ) 

•  cos  (btry) 

(Control  Lat)  =  90  ^  arccos  (x'  1) 

pp  =  great  circle  distance  in  degrees  from  transmitter 
terminal  to  control  point 

4,  GeograpMc  longitude  of  control  points  along  great 
circle  route,  (y  =  y2  *  y') 
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y'  =  arccos 


cos 


«  cos  (90  -  x2) 


Sin  (90 


x3  =  Control  point  geograpkic  latitude 
Si  OeomagnetiG  latitude  of  control  points  along. great 
circle  route,  (y  -  90  -  y' ) 


y'  =  arccos  |"sin  (x4)  •  sin  (x3)  + 
cos  (x4)  •  cos  (x3) 

•  cos  (x3  -  x4)l 


x4  -  latitude  of  geomagnetic  north  pole 

6.  Local  time  at  receiver  terminal. 

LMT  =  GMT  *  y2/l$ 

LMT  e  local  mean  time  hours 
GMT  -  Greenwich  mean  time  -  hours 

7.  Sun' s  zenith  angle  at  control  points. 

cos  (vp)  ^  sin  (z)  .  sin  (ssp)  4  cos  (z)  .  cos  (ssp)  • 
cos  [(15  X  GMT)  -  180]. 

~  sun' s  zenith  angle  <*  degrees 
z  =  control  point  latitude  «  degrees 
ssp  «  latitude  of  subsolar  point  of  sun  for  middle 
of  momh  in  question  «  degrees 
GMT  Greenwich  mean  time  hours 

8.  Ionospheric  absorption  index  "1" 

I  =  (1  +  .  0037  SSN)  (cos  .  881  ^ 
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SSN  4*  predicted  or  observed  12  montb  moving  average 
Zuricb  sunspot  number »  (0  ^  SSN  ^  200) 

9.  E -layer  distance  factor. 

eMc  =  2.08SOOOOOO  •  10‘^  -  1.  31200756  •  10‘^  •  x  + 

3.  761385053  ‘  lO"^  -  5.038476266  ‘  10*^  .  x^  + 

2.624808315  •  lO"'*  x^  -  5.  976618436  •  10 •  x^  + 

1.334494261  •  lO"^  -  4.368460907  .  lO’*^  •  x^ 

X  =  great  circle  distance  -  bundredlUs  of  statute  miles 

10.  £-2000  MUF. 

emc  =  3.345996232  +  37.  67736072  •  1  -  52.41191754  • 

+  39,26151056  •  1^  -  10,66484988  •  1^ 

1  =  absorption  index 

11.  E -layer  MUF. 
emufy  =  (emc)  •  (elfc) 
emc  =  £-2000  MUF 

12.  F -layer  distance  factor. 

fife  =  4.699243101  »  10^^x4  2.26463434)  •  10*^-  x^ 

+  9.202437332  •  lO’^  •  x^  +  6.  865259817  •  I'o"®  •  x"^  - 
9,985831104  •  10’^  .  x^  +  4, 49151441  lO"^  •  yt  - 
6,  712654756  •  lO’^  •  x^ 

X  •  great  circle  distance  -  bundreds  of  statute  miles 

13.  F2-layer  Fourier  generation  of  foF2  and  M-30O0  factors, 
Fourier  time  variation  function,  used  to  obtain  foF2 

or  M-3000  factor: 


ir.'q* 


I  (x.yit)  =  ab  (x, 
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j 


y)  +  (xtV)  cos  (jb  1)  t  + 

jb=2 

ba.,  (x,y)  sia  (jb  =  1)  t  1 


"  I ^is. ka,  io  Sa  ^ 

ka  =  1  jb  -  1,2,3.  i  i .  j 

I 

%  <*•''>  "  I®i,.ka.i6  °ka  ^ 

ka  =  1  jb  =  2, 3, 4. . . .  j 


X  =  geograpbic  latitude  of  coatrol  point  degrees 
(-90*  X  ^  90*) 

y  -  geograpbic  longitude  of  control  point  -  degrees 
(0*  ^  X  ^  360*  *  East  of  Nortb) 

t  =  local  tirne  at  receiver  •  (bour  angle)  (180  ^  t  ^  18Q*) 
I,  j,k,  L  =  Constants  describing  number  of  barmonies 
in  the  Fourier  functions 


io  -  matrix  for  foFZ  or  M-3000  (bigb  or  low  SSN) 

Q.  ,  .  =  Fourier  coefficients  defining  the  function 

is,ka,  lo  - 


(x,y,t)  (Time  variation) 


Latitudinal  and  longitudinal  variation  functions  are 
shown  in  Tables  I  and  2, 

For  a  detailed  explanation  of  the  generation  of  P (x.y.t) 
refer  to  [Jones  and  Gallet,  1961].  The  above  formulas  and 
constants  have  been  altered  from  the  originals  in  tbe  reference 
to  make  tbe  generation  permissible  on  a  larger  variety  of 
computers. 


Table  GeograpMcai  Ft^ctiens  in  Latitude 
and  Longitude. 


MIXED  LATITUDINAL  AND  LONGITUDINAL 


y) 

Cos  (x) 

cos  (y) 

cos  (x) 

sin  (y) 

sin 

(x) 

cos  (x) 

cos  (y) 

sin 

•  • 

(x) 

• 

cos  (x) 

•  ft  • 

sin  (y) 

ft  ft  ft 

•  • 

sin 

• 

La 

ft  ft  ft 

(x)  cos 

ft  ft  ft 

(x)  cos 

(y) 

sin 

La 

(x)  Cos 

(x)  sin 

(y) 

La 

-  k 

-  b 

2 

1 

MIXED  LATITUDINAL  AND  LONGITUDINAL 
VARIATION  -  SECOND  ORDER  IN  LONGITUDE 


cos  (x)  cos  (2  y) 

2 

Cos  (x)  sin  (2  y) 
sin  (x)  cps~  (x)  COS  (2  y) 
sin  (x)  CQS^  (x)  Sin  (2  y) 


sin  (x)  co8“  (x)  cos  (2  y) 
sin  (x)  CQs^  (x)  sin  (2  y) 
Lb  =  I 


14. 
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FZslayer  gyro  frequency  may  be  adequately  represented 

for  MUFbFOT  predictions  by  a  set  of  least  squares 

coefficients  describing  tbe  ortbogonal  polynomial. 

/  ^  ^  2  6,  O' 

=  •••  *1,?*  *••• 

I  ±  ^  2  6,  6 

'^7il  ^7,2*  ^7,3*  •••^7,7*^* 


y  =  gyro  frequency  of  F2  layer  »  Me/  s 
X  -  longitude  of  Control  point  -  degrees 
z  =  latitude  of  control  point  -  degrees 
a  =  Set  of  least  squares  coefficients  describing  the 
ortbogonal  polynomial  (y) 

15.  F2-4000  MUF. 


16. 


17. 


F2-4000  MUF  =  (foF2)  x  {M-4000  factor) 

F‘MUF  for  low  and  bigh  solar  activity. 

F-MUF  =  ZDF  +  fife  •  |r(F2^4000)  -  ZOrj 

ZDF  “  foFZ  +■  I  FZ^layer  gyro  frequency 

M-4000  factor  s  (i,  i)  ^  M-3000  factor 

interpolation  for  intermediate  values  of  solar,  activity. 

(10  £  SSN  ^  180) 

SSN-10)]/170 


F^MUF 


F  -  MUF^^^-  (180 


.  F-MUF^^j( 


F  -  -  F-MUF  for  SSN  10 


F 


F-MUF  for  SSN  180 
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18i  Angle  at  t^e  ionosphere  ((|>). 


<!>'  =  sin  ”^J^P  ^sin  (90  +  tan  ^  (cos  /  2  -  p/  (p  +  i 

~  ^ /2 

<t>'  -  angle  of  incidence  -  degrees 
p  =  radius  of  earth  «  km 
A'  =  great  circle  distance  ■>  degrees 
-  layer  height  «  Ion 

Ionospheric  absorption  (single  reflection) 

A  s  615.  5  N  sec  <}>  (I) 

(f  +  f 

g 

<|>  3  angle  of  incidence  for  100  km  •  degrees 

£  -  operating  frequency 

I  3  absorption  index 

f  =  gyro  frequency  -  Me  /  s 
g 

Basic  transmission  loss  for  isotropic  antennas  in  free  Space. 

Sf  °  *  *°  *'’*10  ^ 

=  36,58*?0Iog,^  H)*Z0  lof!,o(Ve/,> 


p^  =  power  available  at  receiving  antenna 

p  =  power  delivered  to  transmitting  antenna 

d  =  ray  path  distance  «  miles 
\  =  wave  length 

^  frequency  ^  Me  /  s 
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Zli  Relatidnahip  betweeii  9'  aad  A. 
sill  9'  -  GOB  A 


<!>'  =  angle  of  incidenGe  at  ionospliere 
p  =  radius  of  earth 
-  layer  height 

A  =  radiation  angle  of  wave 

Zli  Ground  reflection  factors  for  vertiGal  and  horizontal 
polar  ization^ 


sin  A  *  ^  .  .  2  . 

(«  ^  IX)  «  COB  A 

r 


Sin  A  ^ 


(c  -  ix)  -  cos  A 
r 


(c  •  ix)  sin  A  ^  Vi  .  7  _  . "  .2  * 
if  -  r  (c  »  ix)  =  cos  A 

‘'V  ”  I-:--.-  -:- - - - - - - 


(c  *  ix)  Sin  A  +  VV  i  2  . 
r  (c^  -  ix)  -  cos  A 


c  =  relative  dielectric  constant  of  earth, 

r  ------ 


X  -  18  X  10^  (T/f, 

w  « 

V 

(J  =  conductivity  of  earth  (mhos /meter), 
u)  =  angular  freqviGncy, 
f  =  frequency  in  megacycles, 

A  =  angle  of  elevation  in  degrees*  and 
i  =  V'sl. 

=  dielectric  constant  of  free  space 


! 

■( 


I 


1 


j 


4 


Rbombic  antenna  power  gain  relative  to  isotropic  in  free  space. 


g(A.  0)=  3.  2  ♦  {£)„(A,p>^ 

Oy(A,  0)  =  COS  if  i  ^2  •  sin  1$  •  sin  A 


+  1-2  |K^ I 


cos  (vji^ 


4ek 

X 


sin 


sin  u.  sin  u_ 

Djj  (A,p)  =  cos  (t>  • 

1  2 


sin4> '  GOB  A) 


u. 


r|Kjjl2  ^  1  -  2|Kj^|  cos  ^  •  sin  A)]" 

■j^l  <•  cos  A  •  sin  («j)  +  p) 


cos  A  *  sin  (<|>  p) 


g  (Ai  P)  ^  power  gain  relative  to  isotropic  in  free  space 
A  «  angle  of  departure  in  degrees > 
p  =  angle  of  azinautb  in  degrees, 
i  =  rhombic  leg  len^  in  meters, 
h  =  antenna  height  in  meters, 

<t>  =  tilt  angle  in  degrees. 
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tan 


j  if  •  sin  A  .  Sin  (tan”^*^^  +  •^)  j  + 


p  «  y  .  sin  A 


2  4  2 

|j^  I  ^  [p _ 4  _sin_  A  -  2^_;  sin  _A  _v  jIPiS  c  ] 

^  [p  +  sin^  A  +  2p^  •  sin  A  •  gos  S  ] 


K, 


[2^2  .  4  .  .  .  2  .  -  /  X  7 

ij^p  4  y  .  sin  A  -  2p  .  y  *  sm  A  .  cos  (a  +  2  sin  “i  J 

p  4  y  i  sin^  A  +  2  ^  y^  •  Sin  A  .  cos  4-  sin  ^  ‘^) 

P  ^  2 

y 


*  A  • 

.  Sin  A  •  sm  j 
p  A  Sin  A 


_  2^2 
y  =  €  +  X 

'  r 

ri  2  .  ,2  ^  2  i 

P  '  n  *y  ^  +  X  J  - 


a  =  i-tan  x _ 

(c  »  cos^  A) 
r 

X  =  18  X  10'  cr/f 


=  relative  dielectric  constant  of  tbe  ground 
O'  =  conductivity  of  tbe  ground  in  mhos  per  meter 
f  =  operating  frequency  in  megacycles  per  second 


iKjjl  e^^H 


horizontal  reflection  coefficientt  and 
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Ky  =  “  (Kyi  e^*^V  =  veiftiGal  reflection  coefficient. 


24.  Power  gain  of  liaU«wave  horizontal  dipole. 


=  (.  74582)^  ^co«  ‘  coe  A)  +  1-2  •  K^j  »  cos 

Bin  A 


•'M  K 


•  sin  A)'|^ 


=  amplitude  of  horizontal  reflection  coefficient 

4*0  ~  phase  amplitude  of  horizontal  reflection  coefficient 
n 

H  -  height  of  antenna  -  meters 
\  a  299.  792 5 /frequency  in  megacycles 
Hai 
\ 

A  a  angle  of  elevation 

Note:  72  ohms  assumed  impedance  of  antenna 
®»*''d.cibeU  =  *  *°»10  '®A> 


25.  Power  gain  of  vertical  antennas 


120  r 

"  RA  ’  [• 


cos  (A)  sin  A)  -  cos  (A) 


»  2  ?  K  cos  (v|<„  -  2A  ^  sin  A  | 


=  30  ^  .  |c  +  In  (4A)  «  Cj  (4A)| 


+  I  I  +  cos  (2A)  ]  X  I^C  +  In  (2A)  »  C.  (2A)  | 
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+  tin  (2A)  X  I  (4A)  -  S.  (2A) 

i  2 

G  =  .  S7 72 156649 
H  =  heiglit  of  antenna  -  meters 
Si  =  sine  integral 
Ci  =  cosine  integral 
A  =  2  IT  H/\ 

Efficiency  factor  for  short  vertical  antennas. 

EFF  =  25.  64620146  -  364.8173803  x 

+  2179. 890548  x^  -  6091. 333295  x^ 

+  6416.702573  x'* 

X  =  wave  length  of  anteima 
(1/16  ^  X  ^1/4) 

Oround  reflection  loss. 

2 . 

Relationship  of  field  strength  to  transmission  loss. 
(Decibels  above  isotropic) 

E.  107.2  +  20  1ogjQ(f^^).^ 

f  =  frequency  Me  /  s 
me 

=  transmission  loss  decibels 


START 


Read 

L  Circuit  input  parameters 
2.  idnoBpberic,  geophysieai  and 
Systems  parameters 


X.  generalizeb  block  diagram  of  hf  system  performance  routine 


Calculate 

I.  Great  circle  diStanGe 
Zi  Bearings 

3i  Lpcation  of  control  points 
.4i  LMT  at  Control  points 

5.  Sun' s  zenith  angle  at  control 

points 

6.  Geomagnetic  location  of 

control  points 

7.  Ahsorption  index  I  at  control 

points 

8.  Path  absorption  index  I 

9.  E-layer  control  points  MUFs 
10.  F2-layer  control  points  MUFs 

II.  MUF  and  FOT  for  path 


C  alculate 
F 2 -layer  height 
Distance  per  bop 
Radiation  angle  for  E,  EF2 
and  F 2 -layer  rnodes 
Angle  of  incidence  associated 
with  mode 

Maxinium  E -layer  angle  of 
incidence 

Minimum  F2- layer  angle  of 
incidence 

Dominant  type  of  terrain 
Possible  modes 
Ionospheric  loss  for  path 
Critical  frequency  E-layer 
Critical  frequency  F  2 -layer 
Grotmd  reflection  loss 
Ray -path  attenuation 
Transmitter  anteima  gain 
Receiver  antenna  response 
Mode  with  least  loss 


1.  Set  receiving  antenna 

response  to  isotropic 

2.  Relate  system  loss  to 

field  strength  (dbu) 


Calculate 

1.  LMT  at  receiver  site 

2.  Atmospheric,  man-made 

and  galactic  noise 

3.  Compare  signal  and  noise 


Calculate 

1.  Type  of  path,  e.g. ,  polar, 

short  or  long,  etc. 

2.  Circuit  reliability 


Calculate 

1.  Lowest  frequency  in  Band 
7  with  reliability  of  90% 

or  greater 

2,  Iterate  downward  on  freq, 

until  reliability  is  99% 


Graph  o^  print 
MUF -POT  Option 
Figure  1  or  9 
•■END" 


Print 

System  loss  option 
Figure  5 
"END" 


Field  strength 
Option  Figure  6 


Print 

S/N  ratio  option 
Figure  7 
"END" 


Prmt 

Reliability  option 
Figure  2  or  8 
••END" 


Grai^  or  print 
FOT-LUF  option 
Figure  3,  4  or  10 
■•END" 
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xi.  CONGLUSIONS 

The  computer  solution  of  HF  systems  problems  based  on  estab¬ 
lished  manual  methods  is  efficient  and  practical^ 

The  radio  systems  engineer  may  with  this  basic  tool  check  pre¬ 
dictions  against  operational  data  to  update  idle  prediction  scheme  and 
to  test  new  prediction  parameters  with  little  effort. 

The  routine  is  based  on  average  monthly  values;  therefore,  it  is 
most  useful  in  systems  problems  such  as  allocation  of  fre(^uencies  and 
circuit  design.  Use  of  the  relative  values  produced  by  the  prediction 
scheme  are  no  doubt  more  important  to  the  communication  engineer 
than  are  absolute  values «  The  routine  is  most  valuable  for  the  exper*' 
ienced  engineer  with  an  adequate  knowledge  of  the  shortcomings,  as 
well  as  the  usefulness  of  such  a  prediction  scheme. 
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Instarumenlation.  Engineering  Electronics.  Electron  Devices.  Electronic  Instrumentation.  Mechanical  instru* 
ments.  Basic  Instrumentation. 

Physical  Chemistry.  Thermochemistry.  Surface  Chemistry.  Organic  Chemistry.  Molecular  Spectroscopy.  Ele» 
mentary  Processes.  Mass  Spectrometiy.  Photochemistry  and  Haoiation  Chemistry. 

Office  of  Weights  and  Measures. 

BOULDER,  COLO. 

Ciyogenic  Englneerillg  Lshpratoiy.  Cryogenic  Equipment.  Cryogenic  Processes.  Properties  of  Materials.  CryO': 
genic  Technical  Services. 

CENTRAL  RADIO  PROPAGATION  LABORATORY 

lonomhere  Research  and  Propagatioa.  Low  Frequency  and  Very  Low  Fremency  Research.  lonpaphere  Re^ 
sear^.  Prediction  Services.  Sun-Earth  Relationships.  Field  Engineering.  Hadio  Wwning  Services.  Vertical 
Soundings  Research. 

Radio  Prapafation  Engineering.  Data  Reduction  Instrumentation.  Radio  Noise.  Tropoapheric  Measwements. 
Troposphmic  Analysis.  Propagation-Terrain  Effects.  Radio-Metecrolpgy.  Lower  Atmosphere  Physics. 

Radio  Systems.  Applied  Electromagnetic  Theory.  High  Frequency  and  Very  High  Frequency  Research.  Fre¬ 
quency  lltilization.  Modulation  Research.  Antenna  Research.  Radipdetermination. 

Upper  Atmoaphefe  and  Space  Physics.  Upper  Atmosphere  and  Plasma  Physictk  High  Latitude  Ionosphere 
Physics.  Ionosphere  and  Exosphere  Scatter.  Airglow  and  Aurora.  Ionospheric  Radio  Asnpnomy. 

RADIO  STANDARDS  LABORATORY 

Radio  Physics.  Radio  Broadcast  Service.  Radio  and  Microwave  Materials.  Atomic  Frequency  and  Time-Interval 
Standards.  Radio  Plasma.  Millimeter-Wave  Research. 

Circuit  Standards.  High  Frequency  Electrical  Standards.  High  Fret^ency  Calibration  Services.  High  Frequency 
Impedance  Standards.  Microwave  Calibration  Services.  Microwave  Circuit  Standards.  Low  Frequency  Cdlbratipn 
Services. 


